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^ ; ABSTRACT 

We analyse the Fundamental Plane (FP) relation of 39, 993 early-type galaxies (ETGs) 
Qh' in the optical (griz) and 5, 080 ETGs in the Near-Infrared (YJHK) wavebands, forming 

an optical+NIR sample of 4, 589 galaxies. We focus on the analysis of the FP as a 
£^ ■ function of the environment where galaxies reside. We characterise the environment 

^ ' using the largest group catalogue, based on 3D data, generated from SDSS at low 

redshift (z < 0.1). We find that the intercept "c" of the FP decreases smoothly from 
high to low density regions, implying that galaxies at low density have on average 
lower mass-to-light ratios than their high-density counterparts. The "c" also decreases 
as a function of the mean characteristic mass of the parent galaxy group. However, this 

■ trend is weak and completely accounted for by the variation of "c" with local density. 

The variation of the FP offset is the same in all wavebands, implying that ETGs at low 
density have younger luminosity- weighted ages than cluster galaxies, consistent with 
the expectations of semi-analytical models of galaxy formation. We measure an age 
variation of ~ 0.048 dex (~ 1 1 % ) per decade of local galaxy density. This implies an age 

■ difference of about 32% (~ 3 Gyr) between galaxies in the regions of highest density 
' and the field. We find the metallicity decreasing, at ~ 2 cr, from low to high density. We 

also find 2.5 a evidence that the variation in age per decade of local density augments, 
up to a factor of two, for galaxies residing in massive relative to poor groups. The 
velocity dispersion slope of the FP, "a", tends to decrease with local galaxy density, 
. with galaxies in groups having smaller "a" than those in the field, independent of 

$_i ' the waveband used to measure the structural parameters. Environmental effects (such 

as tidal stripping) may elucidate this result, producing a steeper variation of dark- 
matter fraction and/or non-homology along the ETG's sequence at higher density. In 
the optical, the surface brightness slope, "6", of the FP increases with local galaxy 
density, being larger for group relative to field galaxies. The difference vanishes in 
the NIR, as field galaxies show a small (~ 2.5%), but significant increase of u b" from 
g through K, while group galaxies (particularly those in rich clusters) do not. The 
trend of "6" with the environment results from galaxies residing in more massive 
clusters, since for groups no variation of "6" with local density is detected. A possible 
explanation for these findings is that the variation of stellar population properties with 
mass in ETGs is shallower for galaxies at high density, resulting from tidal stripping 
and quenching of star formation in galaxies falling into the group's potential well. We 
do not detect any dependence of the FP coefficients on the presence of substructures 
in parent galaxy groups. 

Key words: galaxies: fundamental parameters - stellar content - formation - evo- 
lution - galaxies: groups: general 
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1 INTRODUCTION 

A robust prediction of hierarchical theories of galaxy for- 
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in shaping the galaxian properties. From the observational 
viewpoint, this is successfully confirmed by the existence, 
for inst ance, of the w ell established morphology-density re- 
lation (lDresslerl fl98Cn. for which denser environments are 
preferably populated by galaxies with early-type morphol- 
ogy in contrast to the late- type dominated field regions. 

The role of environment in shaping the properties of 
ETGs - the most massive galaxies in the local Universe 
- is still an open question in our understanding of galaxy 
formation and evolution. Studies of the colour-magnitude 
(hereafter CM) relation in the low-redshift-Universe have 
found no significant difference in the average colour of ETGs 
among high- and low-density environments, implying ei- 
ther a tiny difference or a strong anti-correlated variation 
of stellar populati on properties, i.e. age and metallicity , 
with environment i|Bernardi et al.ll2003bl ; iHogg et al]|2004 
lHaines et al.ll200rih . Even small differences in the star for- 
mation history of galaxies residing in different environments 
would be magnified as one approaches the galaxy forma- 
tion epoch. Several studies have found no remarkable dif- 
ference in the properties of cluster a nd fiel d galaxies at 
high redshift. For instance, Koo et al. | (|2005al lbh found that 
both field and cluster E TGs have v e ry sim ilar CM rela- 
tions at redshift z ~ 0.8. ICool et"al] |2006l ) also reported 
that the offset and slope of the CM relation of field and 
cluster galaxies are f ully consistent up to redshift z ~ 0.4. 
IPannella et"al] (|2009h found that ETGs have similar char- 
acteristic ages, independent of the environment they be- 
long to. An opposite picture comes out of spectroscopic 
studies of ETGs at redshift z ~ 0. Many authors have 
reported younger luminosity-weighted ages for field rela- 
tive to cluster galaxies, with the age differ e nce amount- 
ing t o ~ 1-2 Gyr (e.g. iGuzman et alj 19921; Trager et all 



20001; iKuntschner et all 120021; iTerlevich fe Forbes! l2002t 



Thomas et all 120051 ; IBernardi et all 120061 : IClemens et al 



2009). Such difference seems to exist also when con- 



sidering galaxies residing in high-density, low-velocity 

dispersi on systems, suc h as the Hickson Compact 

Groups l|de La Rosa et alll2007l) . 

A even more puzzling issue is the difference in metal 
content of ETGs among different environments, as dif- 
ferent stud ies have found tha t field ETGs are more 
metal-rich (iThomas et all 120051: Ide La Rosa et all 120071 ; 
IKuntschner et al.ll2002l; IClemens et all 1 20091) ? as metal-rich 
as (IBernardi et al. 2006; Annibali et alll2007r ) , or even more 
meta l-poor than their cluster counterparts (|Gallazzi et all 
2006). A main feature of the spectroscopic investigations is 
that the galaxy spectra, and hence the inferred stellar pop- 
ulation properties, refer to the galaxy inner region, typically 
inside one effective radius. On the other hand, ETGs are 
know to possess internal colour gradients, with their central 
part being redder than the out skirts (see IWu et al 1 120051 : 
iRoche. Bernard! fc Hvd"el 120091 ; and refere nces therein) . 
Thes e gradients are mainly due to metallicity (|Peletier et all 
11990^ with a small, but significant contribution from 



age ( La Barbera fc de Carvalh o 2009; Clemens et al.l l2009) 



Moreover, as shown by La Barbera et alH 20051 ). colour gra- 
dients seem to depend on the environment where galaxies 
reside, hence complicating the environmental comparison of 
spectroscopic properties. 

Due to its smal l intrinsic dispersion (see 
e.g. iGargiulo etafl 120091 : iHvde fc Bernardil 120091 : and 



references therein), the Fundamental Plane (FP) relatio n 
of ETGs (|Dressler et all 1 19871 : iDiorgovski fc Davis! Il987h . 
i.e. the scaling law involving radius, velocity dispersion, 
and surface brightness, is a p owerful tool to measure the ir 
mass-to-light ratio (see e.g. Ivan Dokkum fc Franxl Il996l ). 
This results from interpreting the FP as a consequence of 
mass-to-light ratio systematically varying with mass and 
the non- homology of ETGs (see e.g [Hiorth fc MadserJll995l: 
Capelato. de Carvalho. fc Carlberd 1 19951: ICiotti fc Lanzonil 



19971: iGraham fc Collesj Il997l : iBusarello etafl 



Truiillo. Burkert. fc Bell 



Bolton et all 



Tortora et all 120091 ) . Several studies have used 




to constrain the environmental variation, at a given mass, 
of galaxy luminosi ty. Even in thi s case, a contradictory 
picture emerges. iBernardT et all l|2006h found the FP 
relation at low- and high-density to have consistent slopes 
but a significant offset, interpreting it as a difference 
of ~ 1 Gyr in the for mation epoch of field and cluster 
ETGs. On the contrary, iD'Onofrio et all l|2008l ). analysing 
galaxies in massive nearby clusters, found also signif- 
icant variations of the slopes with local density. High 
redshift studies of the FP have reported a consistent 
evolution of the mass-to-light ratio of ETGs among dif- 
ferent environments, implying the same formation epoch 
of field and cluster galaxies 
van Dokk um fc van der Marel 



■fergensen et all 20061: 



van der Wei et all 



(2005) also found that the difference in mass-to-light ratios 
between field and cluster galaxies depends on galaxy mass, 
with low mass s ystems exhibiting a strong differentia l 
evolution (see alsolPahre. Diorgovski. fc de Carvalho! l200ll ; 
Idi Serego Alighieri et alll2006l ). 

This work is the third paper of a series aimed to in- 
vestigate the properties of bright (M r < —20) ETGs as 
a function of the environment where they reside, in the 
low-redshift-Universe (0.05 < z < 0.1). The Spheroid's 
Panchromatic Investigation in Different Environmental Re- 
gions (SPIDER; see La Barbera et al. 2010a, hereafter paper 
I) is based on a sample of 39,993 ETGs with photometry 
and spectroscopy available from SDSS, and Near-Infrared 
(NIR) photome try available from th e UKIDSS-Large Area 
Survey (LAS; L awrence et alj|2007f l. Here, we define the 
environment where the ETGs reside, and investigate the en- 
vironmental dependence of the FP. To this effect, we cross- 
correlate the SPIDER sample to the largest group cata- 
logue, based on 3D data, generated from SDSS at low red- 
shift (z < 0.1). This group catalogue contains 10,124 sys- 
tems selected using a redshift-space friends-of-friends (FoF) 
algorithm from the SDSS data release 7 (DR7, covering 
9,380 square degrees of the northern sky). Notice that we 
chose not to cross-correlate the SPIDER sample with ex- 
isting SDSS-based group catalogues, as they do not cover 
either the entire area or the (low-)redshift range of the 
SPIDER survey. For instance, the maxBCG catalogue is a 
volume limited cl uster sample spann ing the reshift range 
0.10 < z < 0.30 (|Koester et al.ll2007f l. while the SPIDER 
sample is at z < 0.1. The C4 catalogue (|Miller et al.ll2005l) 
spans the redshift range of 0.02 to 0.17, but it covers only 
2, 600 square degrees on the sky, being based on SDSS-DR2. 
The FoF group / cluster catalogue we use in the present study 
has the main advantage of providing a complete sample of 
groups and clusters in the local Universe, covering the large 
sky area of SDSS-DR7. This allows us to derive the FP rela- 
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tion in a wide range of environments using both optical and 
NIR data, and characterising the environments in terms of 
both "local" (e.g. galaxy density) and "global" (e.g. parent 
group mass) observables. The large sample size provided by 
the SDSS+UKIDSS data allows us not only to analyse the 
offset but also the slopes of the FP as a function of environ- 
ment. 

The layout of the paper is the following. In Sec. [2] we de- 
scribe the samples of ETGs. Sec. [3] presents the catalogue of 
galaxy groups we use to define the environment, while Sec. [3] 
describes the properties measured for each group, that de- 
fine the "global" environment of ETGs. Sec. [5] deals with 
the measurement of quantities, such as local galaxy density, 
that define the "local" environment. Sec. [6] introduces the 
FP relation, and the way we measure its slopes and inter- 
cept. In Sec. [7] we derive the intercept of the FP in different 
wavebands as a function of environment. Sec. [8] analyses 
how the offset of the FP constrains the difference in stellar 
population content, i.e. age and metallicity, of ETGs among 
different environments. Sec. [9] describes the dependence of 
the slopes of the FP on local and global environment. In 
Sec. llUI we discuss the main results of this work. A summary 
is provided in Sec. 1111 Throughout the paper, we adopt a 
cosmology with f2 m = 0.3, Q\ = 0.7, and Ho = 75 km s _1 
Mpc" 1 . 



2 SAMPLES OF ETGS 

The SPIDER sample consists of 39, 993 ETGs, with avail- 
able griz photometry and spectroscopy from SDSS-DR6. 
Out of them, 5, 080 galaxies have also photometry avail- 
able in the Y JHK wavebands from UKIDSS-LAS (see pa- 
per I). All galaxies have two estimates of the central ve- 
locity dispersion, one from SDSS-DR6 and an alternative 
estimate obtained by fitting SPSS spectra w ith the soft- 
ware STARLIGHT (|Cid Fernandes et alJliooBT ). In all wave- 
bands, structural parameters - i.e. the effective radius, R e , 
the mean surface brightness within that radius, <fi> e , and 
the Sersic index, n - have been all homogeneously measured 
by the software 2DPHOT (|La Barbera et alj|200gh . 

For the present work, we select two samples of ETGs. 
First, we define an optical, r-band sample of ETGs, con- 
sisting of all galaxies brighter than M r = —20.55, and with 
Xr < 3, where M r is the 2DPHOT total galaxy magnitude 
and Xr is the reduced x 2 of the two-dimensional Sersic fit- 
ting in r-band (see paper I). This selection leads to a volume 
complete sample of ETGs, with better quality structural pa- 
rameters. This optical sample consists of N — 36, 124 ETGs. 
We also define an optical+NIR sample, by selecting only 
the galaxies with photometry available in all wavebands, 
brighter than M r — —20.55, and with x 2 < 3 in all wave- 
bands Q. This selection leads to a volume complete sample of 
4, 589 ETGs, with good quality structural parameters from 
g through K. 

The dependence of the FP relation on the environment 
is first analysed by dividing both samples into field and 



1 Our results are not affected by the choice of the x 2 cutoff value. 
Using a value of 2 would not change at all the environmental 
trends of FP coefficients. 



group galaxies (Sec. l5.2[) . For the subsamples of group galax- 
ies, the role of environment is also analysed with respect to 
group properties, such as mass, galaxy offset from the group 
centre, density (measured relative to other group members), 
as well as substructure (Sec. [5]). As detailed in Sec. 15.21 the 
r-band sample of ETGs (N = 36, 124) is split into two sub- 
samples of 16, 717 group and 11, 824 field galaxies. For the 
optical+NIR sample of ETGs (N = 4, 589), the subsamples 
of group and field galaxies consist of 2, 187 and 1, 359 ETGs, 
respectively. 



3 THE GROUP CATALOGUE 
3.1 The updated FoF catalogue 

The updated FoF group catalogue from SDSS contains 
10, 124 systems and is created as described in Berlind et 
al. (2006), being the difference the area used (9,380 square 
degrees, compared to the original area of 3, 495 square de- 
grees from DR3), as the new version is based on the DR7. 
Besides the coordinates and central redshift, the list gives es- 
timates of the total number of galaxies in the group and the 
velocity dispersion. However, we decided to re-estimate the 
central redshift of the groups, applying the gap technique 
jAdami. Biviano. feMazurd Il9981 ; iLopesi 120071 ; lLopes et ail 
l2009al 'l to the central (0.67 Mpc) galaxies (requiring at least 
three galaxies in this analysis). Doing that we are able to re- 
derive redshifts for 9, 235 (90.4%) systems of the "original" 
sample. Fo r these groups we appl y the "shifting gapp er" 
technique (|Fadda et al.l ll996') as in lLopes et al.l l|2009al ') in 
order to obtain a list of group members independent of the 
FoF algorithm. These groups are then subjec ted to a virial 
analys is analogous to the one des c ribed in | Girardi et al.l 
Jl99St). iPopesso et alJ |2005l , I2007T ), iBiviano et all (|2006h . 
lLopes et al.l ( 2009ah . This procedure yields estimates of op, 
#500, #200, M500 and M200 for most of groups from the 
FoF sample. When performing the virial analysis we require 
the existence of at least five galaxies within the maximum 
radius considered. Out of the 9235 groups we are able to 
measure the above parameters for 8, 083, i.e. 79.2% of the 
initial sample (comprising 10, 124 systems). 

In Fig. [T] we show the output of the "shifting gapper" 
technique. Phase-space diagrams are exhibited for 15 FoF 
groups as an example. The velocity and radial offsets are 
with respect to the group centre. Group members are shown 
with solid squares, while interlopers are plotted with open 
circles. This figure helps illustrating the fact that most sys- 
tems in the FoF group catalogue have low velocity disper- 
sions. 



3.2 Selection Effects: limiting magnitude, 

maximum radius, borders of the survey and 
fiber collision 

To measure properties of the galaxies and the groups where 
they reside, such as velocity dispersion, physical radius and 
mass, care should be taken due to a few selection effects. 
For instance, the limiting magnitude of the survey and the 
maximum radius, R ma x, from which galaxies are selected 
play a crucial role. lLopes et al.l (|2009al ) show that a survey 
complete to M* + 1 is sufficient to derive unbiased estimates 
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Figure 1. Phase-space diagrams of 15 FoF groups shown as examples. The velocity, Vi oa , and radial, R, offsets are with respect to the 
group centre. We apply a shifting gapper procedure for the selection of group members (filled squares) and exclusion of interlopers (open 
circles). Notice that the plots show all galaxies (not only the ETGs) selected from SDSS around each group. 



of velocity dispersion, physical radius and mass. They show 
that missing galaxies in the range M* + 1 < M r < M* + 2 
have no significant impact on cluster parameters. On the 
contrary, missing galaxies at M* < M r < M* + 1 affect the 
estimates of these parameters. For the current work, we only 
use galaxies and groups at z < 0.1. In that regime, SDSS is 
complete to at least M* + 1. 

As discussed by iFadda^ et all (|l996T ) and iGirardi et al.l 



l| 19961 . [ 19981 ^. only values of velocity dispersion obtained 
from all cluster galaxies within a sufficiently large radius can 
be considered as physically meaningful. The choice of Rmax 
may be difficult as one wants to cover the whole cluster, but 
avoid the effects of the large scale structure. This choice may 
also be affected by the area of a pointed obser vation or the 
limits o f a large scale survey. On that matter, lLopes et al.l 
(2009aj) show that a maximum radius of 2.5 h~ l Mpc (3.3 
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Mpc for h = 0.75) is large enough to avoid biased estimates 
of cluster properties. Smaller radius may lead to overestima- 
tion of the velocity dispersion, while larger radius may imply 
a contamination from nearby structures. For large values of 
Rmax, some assumptions done for the virial analysis (see 
below) may not hold. It is also important to stress that the 
magnitude and radial limits will also affect the membership 
selection, and hence the galaxy density estimates (see be- 
low). 

The distance from a group to the edge of the survey also 
has a role when estimating group properties. A very large 
value of Rmax has a large chance of overlapping with a bor- 
der. However, under the assumption of circular symmetry, 
partially missing galaxies in one side of a group do not affect 
the process of selecting members and estimating the veloc- 
ity dispersion. For this work, we set the maximum radius for 
membership selecti on to 7? max = 4.0 Mpc, consistent with 
lLopes et al.l ((2009a). Note that this choice for 7? max and the 
use of clusters at z < 0.1, result in cluster prope rties in good 
agreement to previous resu lts in the literature (|Lopes et al.l 
l2009al V lLopes et all (|2009bh also showed that the scaling re- 
lations derived from these parameters provide excellent mass 
tracers and are in agreement to previously reported results. 

One further issue affects the SDSS spectroscopic sam- 
ple and may have an impact on group parameters, as well 
as galaxy densities. Due to a mechanical restriction, spec- 
troscopic fibers cannot be placed closer than 55 arcsecs on 
the sky. In case of a conflict, the algorithm used for tar- 
get selection randomly chooses which galaxy gets a fiber 
(Strauss et al. 2002). This problem is reduced by spectro- 
scopic plate overlaps, but fiber collisions still lead to a ~6% 
incompleteness in the main galaxy sample. Obviously, this 
issue is more important in high galaxy density regions, such 
as the cores of groups and clusters. Our approach to fix it 
is similar to the one adopted by Berlind et al. (2006). For 
galaxies with r < 18 with no redshifts we assume the red- 
shift of the nearest neighbour on the sky (most of times the 
galaxy it collided with). This may result in some nearby 
galaxies to be placed at high redshift, artificially increas- 
ing their estimated luminosities. Hence, the collided galax- 
ies also assume the magnitudes (in addition to the redshifts) 
of their nearest neighbours, resulting in an unbiased lumi- 
nosity distribution. Notice that this procedure to take the 
fiber collision effect into account might be further improved 
by taking advantage of photo-z and surface brightness in- 
formation of collided galaxies. However, we decided not to 
include this extra information in the correction since (i) the 
fraction of fixed galaxies is quite small (at most 6%, at high- 
est densities); (ii) the above correction procedure is the same 
used to create the FoF group catalogue, and it has shown to 
accurately match the multiplicity fun ction of groups in the 
mock catalogues (|Berlind et al.ll200rj ). In the end, we mea- 
sure densities by fixing the collision effect in a consistent 
way to that used for defining the group catalogue. 

To show that the velocity distribution is not affected 
by the way we apply this correction we show in Fig. [2] the 
velocity offset distribution (relative to group centre) before 
and after fixing the fiber collision issue. The mean values of 
the two distributions, with the corresponding uncertainties, 
are indicated in the bottom left of the figure. 
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Figure 2. Velocity offset (relative to group centre) distribution 
before and after fixing the fiber collision issue. 



4 GROUP PROPERTIES 

In this section we briefly describe the procedures adopted 
for estimating group properties, such as velocity dispersion, 
physical radius (-R200), mas s (M200), and substr ucture. Fur- 
ther details can be found in lLopes et all (|2009ah . These pa- 
rameters are considered when studying the dependence of 
the FP on the environment. 



4.1 Virial Analysis 

For the virial analysis, first we compute the line-of-sight 
velocity dispersion (ap) of all group members within the 
aperture Ra (the radial offset of the most distant group 
member). Ra is normally close to the maximum radius 
(-Rmax = 4.0 Mpc), but could be much smaller in case 
the outermost galaxies are flagged as interlopers. The ro- 
bust velocity dispersion estimate (ap) is given by the gap- 
per or bi-weight estimator, depending on whether the num- 
ber of galaxies available is < 15 (gap per) or > 15 (bi- 
veight; lBeers. Flvnn. fe Gebhardtlll990l) . Follo wing the pre- 



scriptions of Danese. de Zotti. di Tulliol (|l980h . the velocity 
dispersion is corrected for velocity errors. We also obtain an 
estim ate of the projected "virial radius" Rpy l|Girardi et al.l 
1998). A first estimate of th e virial mass is given by (equa- 
tion 5 of iGirardi et"aHll998h 



M v 



ZnopRpv 
2G ' 



(1) 



where G is the gravitational constant and 37r/2 is the de- 
projection factor. 

Next we apply the surface pressure term correction to 
the mass estimate jThe fc Whitdll986h . To this effect, we 
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need an estimate of the R200 radius, which (as a first guess) 
is taken from the definition of Carlberg et al. (1997; equation 
8 of that paper). We assume the percentage errors for the 
mass estimates are the s ame before and aft er the surface 
pressure correction (as in lGirardi et al.lfl998T l. 

After applying this correction, we obtain a refined 
estimate of R200 considering the virial mass density. If 
Mv is the virial mass (after the surface pressure correc- 
tion) in a volume of radius Ra, R200 is then defined as 
R200 = Ra[pv /(200p c (z))] 1 ' /2 ' 4 . In this expression pv = 
3My /(47ri? A ) and p c (z) is the critical density at redshift 
2. The e xponent in this equation is th e one describing 
an NFW (jNavarro. Frenk. fc White3ll997l ) profile near R200 
ijKatgert. Biviano. fc Mazure|[2003 v If we use five-hundred 
instead of two-hundred on this equation, we obtain an es- 
timate of Rsoo- Next, assuming an NFW profile we esti- 
mate M200 (or M500) from the interpolation (most cases) or 
extrapolation of the virial mass My from Ra to R200 (or 
Rsoo)- Then, we use the definition of M200 (or M500) and a 
final estimate of R200 (or R500) is derived. Thi s proce dure 
is analogous to wha t is done by iPopesso et al.l |2007l ) and 
iBiviano et all (120061 ). 

Fig. [3] shows, from top to bottom, the following dis- 
tributions: N200 (number of member galaxies within R200); 
velocity dispersion (in km/s); R200 (in Mpc); and M200 (10 14 
M©). Note that, as expected from the FoF selection func- 
tion, most systems are low richness and low mass systems. 
However, there are still 101 (112) clusters with velocity dis- 
persion (mass) larger than 600 bis" 1 (5 x 1O 14 M0). Six 
clusters have op > lOOOfcm/s, with a maximum value of 
1591 ferns" 1 . 

In a forthcoming contribution of the SPIDER series, we 
will analyse the scalin g relations of grou ps and clusters in 
the FoF catalogue (cf. lLopes et al.ll2009bl ). 



4.2 Substructure estimates 

On what regards substructure, we employ a 3D substructure 
test to all FoF groups with at least five galaxy member s 
within R200. Th e PS, or A test jDressler fc Shectmanlll98Sl ). 
is applied as in lLopes et al.l (|2009a ). The code used for th e 
substructure analysis is the one from iPinknev et all l|l996l) . 
who evaluated the performance of thirty-one statistical tests. 

The algorithm computes the mean velocity and stan- 
dard deviation (a) of each galaxy and its N nn nearest neigh- 
bours, where N nn = N 1//2 and N is the number of galaxies 
in the group region. Then these local mean and a values are 
compared with the global mean and a (based on all galax- 
ies). For every galaxy, a deviation from the global value is 
defined by the formula below. Substructure is estimated with 
the cumulative deviation A (denned by ^ Si; see below). For 
objects with no substructure we have A ~ N. 



j-2 / N nn + 1 \ <(- -n2 , , si 

Oi = I ] [(Vlocal —V)+ (aiocal - O) 



(2) 



Any substructure test statistic has little meaning if not 
properly normalised, which can be achieved by comparing 
the results for the input data to those for substructure-free 
samples (the null hypothesis). For the A test the null hy- 
pothesis files are created after the velocities are shuffled 



randomly with respe ct to the positions, which remain fixed 
jPinknev et al.lll996l ). 

For each input data set, we generate five-hundred sim- 
ulated realisations. We then calculate the number of Monte 
Carlo simulations which show more substructure than the 
real data. Finally, this number is divided by the number of 
Monte Carlo simulations. We set our significance threshold 
at 5%, meaning that only twenty- five simulated data sets can 
have substructure statistics higher than the observations to 
consider a substr ucture estimate signif ic ant. Further details 
can be found in IPinknev et al] || 19961 ). lLopes et al l (|2006l ) 
and lLopes et al.l (|2009al ). 

For each FoF group the substructure test is applied to 
the 3D distribution of all galaxy members, as selected by the 
shifting-gapper method. The A test is only applied to sys- 
tems with at least five galaxies available within R2oo- Out 
of the 8,083 groups satisfying these criteria, we find that 
3,570 (~ 44,2%) show significant signs of subs tructure. We 
note th at this fraction is larger than found by lLopes et al.l 
|2009al ). That is probably due to the fact that (in order to 
correct for the fiber collision issue) we artificially assign red- 
shifts to galaxies not targeted for spectroscopic observations. 

We found that the fraction of groups with substructures 
increases from ~ 30% at A/200 ~ 10 1 3M a to ~ 80% for 
clusters as massive as 1O 14 M0. Notice that this result should 
be taken with caution, as low (relative to high) mass systems 
usually have a lower number of member galaxies and this 
might affect the efficiency in detecting substructures at the 
low mass regime 0. 



5 THE ENVIRONMENT OF ETGS 
5.1 Density Estimates 

Galaxy density is estimated as follows. For each galaxy mem- 
ber of a given FoF group, we compute its projected distance, 

1/2 

djv, to the Nth galaxy of that group, where N = N ' mi) , and 
N mem 6 is the number of member galaxies in the group. The 
local density Ejv is given by N/nd%, measured in units of 
galaxies/Mpc 2 . The djv is expressed in units of Mpc, using 
the redshift of the given group. For the density estimation 
we only consider member galaxies within a fixed luminosity 
range (see comments below). Notice that differently from 
other works density is not computed relative to all other 
galaxies in the field, but only relative to the other group 
members. Hence, instead of using £5 or £10, which are more 
common in the literature, we use £jv so that the distance 
used for estimating density scales according to the number 
of members (or group mass). A further note is that if for 
a given galaxy there are fewer than N member neighbours 
(in the given luminosity range), we cannot compute the dis- 
tance to the Nth galaxy of that group, djv. In that case, £jv 
is set to a lower limit of 0.01. However, that only happens 
for about 0.1% of the objects in the optical sample of ETGs 

with gr oup membership. 

In lLa Barbera et ah! (|2010bl) (hereafter paper II), we 
showed that the sample of ETGs is complete to M r = 



2 This issue will be addressed in the forthcoming contribution 
where we analyse the scaling relations of groups and clusters in 
the FoF catalogue. 
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—20.55 for 2DPHOT total magnitudes. This corresponds to 
a limit of M r = —20.32 for SDSS model magnitudes. Here 
we compute local density by (SDSS) petrosian magnitudes. 
The above limit translates to M r — —20.23 for petrosian 
magnitudes. The density estimates are derived in the same 
luminosity range for all galaxies, no matter which FoF group 
they belong to. This limit is fixed by the maximum red- 
shift of SPIDER (z = 0.095). At this redshift, for petrosian 
magnitudes, the characteristic magnitude of the galaxy lu- 



minosity function is M* — —21.43 in r band (Popesso et 
al. 2005). If we make M* + X = -20.23, we find X = 1.2. 
Hence, when deriving densities, we select only group mem- 
bers brighter than M* + 1.2. 

Fig. [4] shows the comparison of Ejv computed with the 
original spectroscopic SDSS survey and after fixing it for the 
fiber collision issue (see section 2.2). In this Figure, we use 
dots to show the results for all galaxies for which density 
was estimated. Large green points show the median values 




Figure 4. Comparison of Sjy computed with the original spectro- 
scopic SDSS survey and after fixing it for the fiber collision issue. 
Dots are used for all the galaxies for which density is estimated. 
Large green points show the median values for bins containing 
5000 galaxies each. 

for bins containing 5, 000 galaxies. As expected, the correc- 
tion increases the Ejv value. The effect is slightly stronger 
at high E m , the regime where we expected the fiber collision 
issue to be more relevant. Fig. [S] shows the density distribu- 
tion for all member galaxies within the 8, 083 FoF groups. 
Galaxies residing in the FoF groups span a wide range of 
environments, corresponding to more than three decades in 
local galaxy density. 

5.2 Definition of group and field galaxies 

The environment of ETGs is first defined by the location of 
a galaxy in either a group or the field. For galaxies residing 
in groups, the environment is further characterized in terms 
of the group properties, such as the group-centric distance, 
mass, physical radius, and substructure. Group galaxies have 
also local density estimated relative to other group members 
(as described above), while for field galaxies, we use the FP 
relation itself to introduce a fictitious local density value (see 
Sec. ED - 

Galaxies residing in groups are simply those having 
group membership according to the shifting gapper method 
(within ii m ax = 4.0 Mpc; see Sec. [3|. Both the optical and 
optical+NIR samples of ETGs (Sec. are cross-correlated 
to the list of galaxies belonging to groups, so that we define 
which SPIDER ETG has a parent group in the FoF updated 
catalogue. 

Field ETGs are defined as follows. We start selecting 
only the ETGs that are not members of any group in the 
updated FoF catalogue. As reported in Sec. [3] because of 
the low number of galaxies available for some (poor) groups, 



Figure 5. Distribution of local galaxy density (£;v) for all group 
members. 

we are able to measure group properties and define galaxy 
members for 79.2% of the FoF groups. This implies that 
some galaxies with no assigned membership may still reside 
in one of the not-measured groups, and should be excluded 
from a genuine field sample. To clean the field sample from 
these cases, we check all the FoF groups (10,124 systems), 
removing those ETGs with radial offsets < 5 x R pe rp,rms 
and velocity offsets c x A z < 5 x a v from any group, where 
Rper P ,rms and er„ are the projected rms radius of the group 
and v elocity dispersion a s defined in the original FoF cata- 
logue l|Berlind et al . 2006). Galaxies on the border of a group 
in the phase-space diagrams (Fig. [1} can be not classified as 
group members by the shifting gapper technique, but never- 
theless their properties might still be affected by the group 
because of the in-falling in its potential well. To remove this 
further source of "contamination" from the field sample, 
we check once more all the 8,083 measured groups, remov- 
ing objects with radial offsets < 5 x R200 and velocity offsets 
c x A z < 5 x crp from any group, where -R200 and op are 
the characteristic radius and velocity dispersion of a group 
computed as described in Sec. [4] The remaining ETGs define 
the field sample. 

Applying this procedure to the r-band sample of ETGs 
(N — 36, 124), we obtain two subsamples of 16,717 group 
and 11, 824 field galaxies. For the optical+NIR sample of 
ETGs (N = 4,589), the above selection results in 2,187 
group and 1,359 field galaxies, respectively. 



6 MEASURING THE FP RELATION 

As in paper II, we write the FP relation as: 

logiie = a log (To + b </Jb> e +c, (3) 
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where "a" and "b" are the slopes, and "c" is the offset (in- 
tercept) of the FP, log R c is the logarithmic effective radius 
of galaxies, in units of kpc, < /i > c is the mean surface 
brightness within that radius, and <to is the central veloc- 
ity dispersion. For the present study, we use ao values from 
SDSS, although we verified that results are the same when 
using STARLIGHT velocity dispersions (see paper I). The 
cto's are corrected to a circularised aperture of r adius r e /8, 
following I j0rgensen. Franx. fc Kiaergaardl (|l995 ). The com- 
putation of log Tie and </i> c in the different wavebands is 
described in detail in paper I. In general, we refer to "a" , 
"b", and "c", as the coefficients of the FP. 

In order to analyse environmental effects on FP coeffi- 
cients, we proceed as follows. In Sec. [7] we first compare the 
offset of the FP in different environmental bins. The values 
of "c" for samples of ETGs in the different bins are com- 
puted by fixing the slopes of the FP to the values of a — 1.39 
and b — 0.315 obtained in paper II for the entire SPIDER 
sample. This allows us to interpret differences in "c" as dif- 
ferences in the average mass-to-light (hereafter M/L) ratio 
of ETGs among different environmental bins (Sec. [8]). No- 
tice that this FP-based M/L ratio differs from the quan- 
tity M dyn /L, where dynamical mass0, M dyn oc Re&o, and 
galaxy l uminosity, L, are dire ctly estimated from the data 
(see e.g. iBernardi et a"i1l200 3a'). The "c" is measured as 



c = log Re — a X log (To — 6 X < fi > e , (4) 



where log ii e , log <n> e , and logtro are the median values 
of logi? e , </!> c , and logcro- For a given sample of ETGs 
(i.e. a given environmental bin), the uncertainties on "c" 
are computed by the width of the distribution of "c" values, 
as computed by Eq. 2] from 1000 bootstrap realisations of 
the sample. Since the values of "a" and "b" are kept fixed in 
Eq.[4]for all realisations, this procedure neglects the effect of 
the (correlated) uncertainties on "a" and "6". As discussed 
in Sec. this further source of uncertainty can be completely 
neglected, as far as the difference of "c" among different 
environments is concerned. The environmental variation of 
"c" is used in Sec. [9] to properly measure the slopes of the 
FP for the ETG's samples populating the different (envi- 
ronmental) bins. The values of "a" and "b" are estimated 
by an orthogonal robust fitting procedure and corrected for 
selection effects as well as the effect of correlated uncertain- 
ties on structural parameters (see paper II for details). We 
implement a procedure to account for the different average 
mass-to-light of galaxies in the different environmental bins, 
as well as the "luminosity segregation" of galaxy populations 
with the environment (see Sec. [9] for details). The scatter of 
ETGs around the FP and its dependence on environment is 
not analysed here, being the subject of a forthcoming paper 
in this series. 



3 In order to determine the proportionality factor relating M dyn 
and -R e a-Q, one has to assume a given galaxy model, making some 
assumptions about the (three-dimensional) distribution of dark- 
and luminous matter in a galaxy. On the other hand, the FP- 
based M/L does not have to rely on any of these assumptions. 
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Figure 6. Variation of the offset "c" of the FP as a function 
of local galaxy density, Ejv. Filled circles plot the values of "c" 
obtained binning the group sample of ETGs with respect to Sjv, 
adopting 30 bins with the same number of galaxies each. Error 
bars denote 1 cr uncertainties. The "c" has been computed by fix- 
ing "a" and "6" in all bins to the same values of 1.39 and 0.315, 
respectively, as shown in the upper-right corner of the figure. The 
dashed black line is the best-fit line to the filled circles, derived 
by a least-squares fitting procedure with "c" as dependent vari- 
able. The slope value of the line is reported in the lower-left of 
the plot. The value of "c" for the field sample of ETGs is also 
shown in the figure as an empty circle. The average local density 
value of field galaxies, marked by the vertical dotted line, is as- 
signed by using the above linear fit (see the text). Red and blue 
squares are obtained by selecting only galaxies in a narrower red- 
shift range and in the lower contamination sample of ETGs. The 
corresponding least-squares fits of "c" versus log E jv are plotted 
by the red and blue dashed lines, whose slope values are reported 
in the lower-left of the Figure. 

7 THE OFFSET OF THE FP AS A FUNCTION 
OF ENVIRONMENT 

We start analysing the dependence of the offset "c" of the FP 
on galaxy environment. Sec. l7.1| presents the results obtained 
from the optical sample of ETGs, while the variation of "c" 
from g through K is analysed in Sec. 17.21 

7.1 The FP offset in r band 

Fig.|6]plots the value of "c" for the optical sample of ETGs, 
as a function of the local galaxy density Ejv. The sample of 
group galaxies is binned with respect to Ejv, with each bin 
including the same number of 557 galaxies. This leads to 30 
density bins. The Figure shows that "c" smoothly decreases 
from low to high density regions. A linear fit of this trend, 
with 

c = a Ejv + fcsjv log Eat , (5) 

gives a Sjv = -8.8734 ± 0.0018 and 6 Sjv = -0.014 ± 0.002. 
The value of b^ N differs from zero at ~ 7 a, meaning that 
the variation of "c" with Ejv is highly significant. The quan- 
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Table 1 . Covariancc matrix terms of the uncertainties 
on a, b, a Ejv , and b Ejv . 



X 


Y 


C{X,Y) 


a 


a 


4.150e - 04 


a 


b 


9.964e - 06 


a 


a Sjv 


-1.515e - 03 


a 




-2.900e - 06 


b 


b 


1.970e - 06 


b 




-7.462e - 05 


b 


bs N 


-3.499e - 08 


a Sjv 


«s» 


6.378e - 03 


a Sjv 




7.565e - 06 






5.137e - 06 



tity 6ejv represents the variation of "c" per decade of local 
galaxy density, and encodes all the information about the 
dependence of the FP offset on the (local) environment (see 
Sec. [SJ. For the field sample of ETGs, an estimate of the 
local density is not available. Hence, we assign a fictitious 
Ejv value, "Efuid, to all the field ETGs. To this aim, we first 
compute the "c" for the field sample. Then, inverting Eq. 
with as N and bs N given by the values listed above, we infer 
the S/ieid value. 

A possible source of concern in the trend of "c" with en- 
vironment is the contamination of the SPIDER samples from 
galaxies with non-genuine early-type morph ology (see pape r 
I). Due to the morphology-density relation (|Dresslerlll98d) . 
these objects are expected to be more abundant in the low 
density regions. To analyse the impact of such contamina- 
tion, in paper I we have defined a low contamination sample 
of ETGs consisting of 32, 650 (out of 39, 993) galaxies. For 
each bin of local density, we select only galaxies in this low 
contamination subsample. The resulting trend of "c" with 
local density is shown by the blue symbols in the Figure. Be- 
sides a rigid shift, the slope of the trend is identical to that 
obtained with the entire sample, proving that contamination 
does not affect at all our results. Another source of concern 
is the redshift range of the SPIDER sample. Fainter galaxies 
are relatively more abundant in the lower density regions, 
and are better represented in the low redshift tail of the sam- 
ple. Due to small evolutionary ef fects in the luminosity of 
ETGs (e.g. iBernardi et al.ll2003ar ). this effects can bias the 
trend of the offset with the environment. To tackle with this 
issue, we select in all the environmental bins only galaxies 
in a narrower redshift bin, from z = 0.0625 to z = 0.0825. 
The corresponding trend of "c" with local density is shown 
by the red symbols in Fig. [6] The slope of the "c"-logEjv 
fit is fully consistent with that of the entire sample. 

We find that the "c"-logEjv trend is robust with re- 
spect to the statistical estimator used to calculate the 
FP intercept. Replacing the median values of logi? c , < 
u > e and log <7o, in Eq. U with the peak values of 
the distributions, as estimated by the bi-weight estima- 
tor l|Beers. Flvnn. fc Gebhardtlll990D . we obtain a very sim- 
ilar trend, whose slope, &e n , is indistinguishable from that 
reported in Fig. [6] Since the bi-weight (relative to the me- 
dian) estimator is less sensitive to the shape of the tails in 
the distribution being analysed, this test also shows that 
the trend of "c" with log Eat is not affected by the different 
shape of the massive end of the galaxy luminosity function 
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Figure 7. Same as Fig. \E\ but plotting "c" as a function of the 
normalised cluster-centric distance R/ -R200 • The trend of increas- 
ing "c" with R/R200 closely reflects the one shown in Fig. [6] for 
the local galaxy density. Notice that the field sample is also plot- 
ted in the Figure as an empty circle. The cluster-centric distance, 
R f ield/ R200, is assigned to field galaxies in the same way as the 
local galaxy density, Ey^^ (see the text). 



in different environments, i.e. the larger fraction of bright 
galaxies inhabiting higher (relative to lower) density envi- 
ronments. Notice that the values of "c" in different bins of 
log Ejv are measured by keeping the values of FP slopes, 
"a" and "6", fixed (Sec. [6}. On the other hand, as shown in 
Sec. [9] the FP slopes show a small, but significant trend with 
the environment, which might bias the "c"-logEjv trend. 
In order to address the issue, we recomputed the "c" by 
setting a — 2 and b — 0.4 in Eq. [3] This corresponds to 
measure the mass-to-light ratio of ETGs by means of the 
virial theorem, rather than the FP equation, hence avoid- 
ing the possible "degeneracy" between the environmental 
variation of FP slopes and intercept. Setting a = 2 and 
b = 0.4, the trend of "c" with log Ejv turns out to be fully 
consistent with that shown in Fig. [6] with the slope of the 
best-fitting line amounting to —0.02 ± 0.004 (compared to 
—0.014 ± 0.002). More in general, we verified that all the 
environmental trends of "c" remain statistically unchanged 
when using either the FP or the virial theorem equation. 

The error bars in Fig. [6] deserve further comment. As 
stated in Sec. [5] we compute the uncertainties on "c" by 
keeping "a" and "b" fixed in Eq. [4] This corresponds to 
neglect the effect of (correlated) errors of "a" and "b" on 
the error bars of "c" in Fig. [5] and hence on the estimated 
uncertainty on bs N ■ To analyse the effect of the errors on 
"a" and "b" , we adopt the following approach. We take the 
covariance matrix of the uncertainties on "a" and "6", as es- 
timated by the orthogonal fitting procedure (see paper II), 
and perform 1000 iterations, shifting each time the values of 
"a" and "b" according to their correlated errors. For each 
iteration, the values of "c" in all the environmental bins of 
Fig0are re-computed, together with the values of oej, and 
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bs N . Notice that in a given bin of Eat, the "c" is computed 
by extracting a bootstrap realisation of the galaxy sample 
in that bin. We then estimate the covariance matrix of the 
values of "a" , "b" , as N , and bs N , among the different itera- 
tions. The uncertainties provided by this covariance matrix 
account for both the scatter around the FP, as already given 
by the error bars in Fig. [7] as well as the correlated errors 
on FP slopes. All terms of the covariance matrix, C(X, Y), 



where X (Y) is one of the quantities 



/ r // 

, 



and 



bs N , are reported in Tab. [T] The root square of the terms 
C(a,a) and C(b,b) are the uncertainties on "a" and "b" 
in r-band (cf. tab. 6 of paper II). The 1 a uncertainty on 
fesjv amounts to ~ 0.0023, being given by the root square 
of C(6s JV , &Sjv )(= 5.137e — 06). This uncertainty is indistin- 
guishable from the value of 0.0022 reported in Fig. where 
the uncertainty on bs N was estimated by neglecting the un- 
certainties on "a" and "b". In other words, as far as the 
environmental trend of "c" is concerned (i.e. the value of 
bs N ), the effect of the errors on "a" and "b" can be com- 
pletely neglected, and the dominant source of error is the 
rms of residuals around the FP. Notice also that the covari- 
ance term between bs N and "a" is negative, implying that 
as "a" increases the value of bs N tends to become smaller. 
The same holds for the covariance term C(b,bs N )- This is 
consistent with the fact that, when fixing the values of "a" 
and "b" to the virial theorem expectations, we obtain a value 
of 6ejv slightly smaller, but still statistically consistent (see 
above), than that reported in Fig. 0. 

Fig. plots the "c" in r band as a function of the me- 
dian cluster-centric distance of ETGs to their parent groups. 
Cluster-centric distances are normalised to the correspond- 
ing -R200 values, allowing a meaningful comparison of galax- 
ies residing in groups having different physical radii. The 
"c" significantly increases from about —8.9 in the very cen- 
tral cluster regions, R/R200 ~ 0.1, to about —8.865 in the 
group's outskirts (R/R200 ~ 7). The range of variation is 
very similar to that found for local density, implying that 
both trends are providing the same kind of information. It 
is interesting to notice that the linear fits of "c" vs. logEjv 
and R/R200 have very similar rms values, amounting to 
0.009 ± 0.001 and 0.008 ± 0.001, respectively. Notice that 
Fig. also plots the value of "c" for the field sample. To 
this aim, a fictitious cluster-centric distance, Rfieid/ R200, is 
assigned to field galaxies, using the linear fit of "c" versus 
R/R200, in the same way as for E^^ (see above). 

Fig. [S] plots the "c" as a function of the median mass 
value, M group, of the groups where galaxies reside. We find 
that the offset of the FP decreases as the parent halo mass 
increases. The variation is significant at the 3 a level, as 



shown by the linear fit of "c versus M g 



whose slope 



value, 0.053 ± 0.018, is reported in the same Figure. The 
variation of "c" with group mass is much weaker than those 
found for R/R200 and Eat. This is also shown in Fig. [9l 
where we plot the "c" as a function of the projected mass 
density, M group / 'R200, of parent groups. This plot can be 
more directly compared to the variation with local density 
Ejv. The FP offset tends to decrease as the "global" density 
increases. However, in the latter case, the variation is much 
weaker than that with Eat. Using Eq. we can correct the 
values of "c" in the different bins of M group / 'R200, removing 
the expected variation due to the different median value of 
Ejv of galaxies in them. The corrected points are shown as 
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Figure 8. Same as Fig. [6] but plotting "c" as a function of the 
mass of the parent group where galaxies reside. The field sample 
is not shown in the plot. 
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Figure 9. Same as Fig. [6] but plotting "c" as a function of 
the projected mass density, M group / -R^ou' °^ *he P aren '- groups 
where ETGs reside. Red circles are the values of "c" corrected 
for the expected amount of variation of "c" among the different 
bins due to the different local density of galaxies in them. The 
corresponding best-fit line is shown by the red dashed line. Its 
slope value is reported in the lower-left of the plot. 



red circles in the Figure. The trend of "c" with "global" mass 
density disappears, proving that the environmental variation 
of the FP intercept is essentially driven by Ejv. In other 
terms, local environment turns out to be the main driver of 
differences in FP offset. 

In order to further analyse this point, we derive the 
variation of "c" with density for different bins of M gr0 u P - 
The r-band sample of ETGs is binned simultaneously with 



respect to Ejv and M g 



with a total of 64 bins, each bin 
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Figure 11. Slope of the "c" vs. log Eat relation. as a func- 

tion of the median mass of parent groups where ETGs reside. The 
slope of the best-fit line (dashed line in the plot) is reported in 
the lower-left part of the plot. The slope value differs from zero at 
the 2.5 a significance level. The dot-dashed lines mark the ±1 a 
confidence levels around the best-fit. 



including the same number of 261 galaxies. Fig. [T0]plots the 
"c" as a function of Eat for different bins of group mass. 
The "c" tends to decrease as Ejv increases in all mass bins. 
Each trend is modelled by a linear fit as in Fig. [6] (see Eq.[5]), 
and compared to the linear fit obtained for the entire sample 
(dot-dashed lines in the Figure). The slope value, bs N , of the 
"c" vs. log £jv relation tends to become progressively steeper 
from poor groups to rich clusters. In fact, the bs N value 
changes from about —0.005 for M group ~ 0.1 x 1O 14 M to 
about —0.023 for M group ~ 5.8 x 1O 14 M0. For comparison, 
the slope of the "c" - log Ejv relation for the entire r-band 
sample is b^ N = -0.014 ± 0.002 (see Fig. EJ). Fig, fill plots 
the &Ejv as a function of the median group mass M group . A 
simple linear fit shows that the steepening of "c" - log E at 
relation with group mass is significant at the 2.5 a level. 



7.2 The FP offset from g through K 

Using the optical+NIR sample of ETGs, we derive the vari- 
ation of the FP offset as a function of environment, from 
the g through the K band. We bin the optical+NIR sam- 
ple of group galaxies with respect to Eat, each bin includ- 
ing the same of number of 273 galaxies. This results into 
a total of eight bins. For a given bin, the "c" is computed 
from Eq. 2] using the same sample of galaxies in all wave- 
bands. Fig. [l2]plots the FP offset as a function of Eat from g 
through K. The values of "c" for the optical+NIR sample of 
field galaxies are also shown in the Figure, adopting the Eat 
value (E f^id) derived for field galaxies in Sec. 17.11 The "c" 
smoothly decreases from low to high density regions, with 
all wavebands showing the same trend. In order to charac- 
terise the trends, we fit the data in each band with a linear 
relation: 

c = a s N ,x + fe Ejv ,x lo S E iv. (6) 



Table 2. Coefficients of the linear fits to the trend of 
"c" with local density in different wavebands. 



waveband 






9 


-9.1452 ± 0.0030 


-0.0124 + 0.0034 


r 


-8.8713 ± 0.0024 


-0.0097 ± 0.0027 


i 


-8.8244 ± 0.0025 


-0.0096 ± 0.0027 


z 


-8.7134 + 0.0026 


-0.0108 + 0.0025 


Y 


-8.5496 ± 0.0027 


-0.0098 ± 0.0028 


J 


-8.5357 + 0.0026 


-0.0080 ± 0.0028 


H 


-8.3709 ± 0.0027 


-0.0119 + 0.0030 


K 


-8.2286 ± 0.0027 


-0.0153 + 0.0030 



with X = grizYJHK . The fit is done by a least-squares 
method, with "c" as dependent variable. The uncertainties 
on fitting coefficients are estimated by TV = 500 bootstrap 
iterations, shifting the values of "c" according to their un- 
certainties. The best-fit values of a», v and „ are re- 
ported in Tab. [2] Notice that the variation of "c" per decade 
in local density, i.e. the b s ^ x , are fully consistent, within 
the corresponding uncertainties, among the available wave- 
bands. The values of a Ejv r (= -8.8713 + 0.0024) and 6 Ejv r 
(= -0.0097 + 0.0027) are fully consistent with those derived 
for the optical sample of ETGs, i.e. a Ejv = -8.8734 + 0.0018 
and 6 Ejv = -0.014 + 0.002 (see Sec.[7T]). We do not attempt 
here to derive the "c" as a function of both Eat and M group . 
In fact, due to the smaller sample size of the optical+NIR 
sample of ETGs, binning the "c" with respect to both Eat 
and M group would overly reduce the sample size in each bin. 



8 STELLAR POPULATION PROPERTIES OF 
ETGS IN DIFFERENT ENVIRONMENTS 

Using the virial theorem, one can rewrite the equation of the 
FP as a power-law relation between the mass-to-light ratio 
of ETGs and two other variables, such as the mass M, lu- 
minosity L, (To, R e , and mean surface brightness (see paper 
II). The intercept of the FP turns out to be proportional 
to the mean M/L ratio of ETGs at a given point of the 
galaxy sequence, parametrised in terms of the two adopted 
variables. Hence, the strong variation of "c" with local envi- 
ronment, (Sec. l7.ll and Sec. l7.2[l . can be actually interpreted 
as a variation of M/L with environment. Such a variation 
can be due to a change of M, L, or both quantities with Eat. 
In order to analyse this point, we plot in Fig. [13] the peak 
value of the quantity \ogR e — alogoo, i-e- the combination 
of log R a and log cro that enters the FP, as a function of 
local galaxy density, for the r-band sample of ETGs. The 
peak value is computed by using the bi-weight estimator. 
We see no significant trend with respect to Eat, implying 
that the trend of "c" is due to a variation of mean surface 
brightness with local density. Fig. Q3] also plots the quantity 
log Re + 2 log <To , which is a proxy for the galaxy dynamical 
mass, M, as a function of Eat. Even in this case, no vari- 
ation with the environment is detected. This suggests that 
the "c"-log Eat relation origins from a change in the aver- 
age luminosity, i.e. stellar population content, of ETGs as a 
function of environment, rather than a variation of galaxy 
mass (see also Sec. I10p . The optical and NIR data can help 
us to further constrain the origin of such variation in stel- 
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Figure 10. Offset "c" of the FP as a function of local galaxy density for different bins of parent group mass, M g 



Each panel is the 



same as Fig. |U but plotting only the FP offset for ETGs in a given bin of M group . The median value of M group , <M gr 

oup^i increases 

from left to right and top to bottom, being reported in the lower-left corner of each panel. The dashed lines are the linear fits of "c" 
vs. logSjv, with their slope values being reported in the upper part of each panel. The dot-dashed line is the same for all panels, and 
corresponds to the best-fitting line of "c" vs. log Sjv for the entire r-band sample. It coincides with the dashed black line in Fig. [6] 



lar content. Assuming for simplicity that the variation of 
M/L with environment can be entirely ascribed to age and 
metallicity, we write the following equations: 

Sc x _ SlogM/Lx 
~~ <51ogEjv ~ (HogEjv 
8 log M*/L x <5 log* 
8 log t 8 log E jv 

6\ogM*/L x SlogZ 

SlogZ 51ogEjv' 1 ' 

where M/Lx (M+/Lx) is the dynamical (stellar) mass- 
to-light ratio in the waveband X (= grizYJHK), 
the quantities ct t x = 8 log M t /Lx/S log t and cz,x = 



5 log M+/Lx /Slog Z are the partial derivatives of M*/Lx 
with respect to the age t and metallicity Z. Since 
SlogM/Lx = 8 log M/M* + 8 log M^/Lx, Eq. [7] assumes 
that, on average, the quantity log M/M*, i.e. the aver- 
age fraction of dark- to stellar-matter in ETGs, does not 
change significantly with the environment. As discussed 
in Sec. 1101 this assumption is supported by the fact that 
the typical dynamical-mass of the ETGs in our sam- 
ple does not change significantly with environment (see 
Fig. 1 14 p . and a variation of M+, at fixed M, is unlikely.. 
The quantities Slogt/Slog Ejv and 8 log Z/5 log Ejv repre- 
sent the variation of age and metallicity per decade in lo- 
cal galaxy density, and are the parameters we aim to con- 
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Figure 12. Dependence of the FP offset on local galaxy density 
in the grizY J H K wavebands. Different wavebands are plotted 
with different colours, as shown by the labels below the dashed 
lines in the left part of the plot. Galaxies in the optical+NIR 
samples of ETGs residing in groups are binned according to Ejvj 
with each bin including the same number of galaxies. For each 
bin, the "c" is derived for the same sample of galaxies from g 
through K. The values of "c" for the field sample (empty circles) 
are reported in the same figure, adopting the local density value 
assigned to field galaxies from the linear fit in Fig. [6] (see Sec. 17. 111 . 
The dashed lines are the linear fits to the trend of "c" with Sjv 
in each band. Error bars denote 1 a uncertainties. Notice that 
suitable shifts have been applied to the values of "c" in different 
bands to better display all the trends in the same plot. 



st rain. We estimate the quantities C t.x and cz.x as detailed 
in lLa Barbera fc de Carvalhol (|2009h and paper II. In short, 
we const ruct several simple stellar population (SSP) models 
with the iBruzual fc CharlotJ l)2005t ) synthesis code. All mod- 
els have a Scalo IMF, and different ages and metallicities. 
Age values span a range of ~ 8 to ~ 13 Gyr, while metallicity 
varies from 0.2 to 2.5Zq. For each model (i.e. each age and 
metallicity), we compute the mass-to-light ratios, M+/Lx, 
by folding the corresponding SED with the grizYJHK fil- 
ter curves. For each band, the values of M+/Lx from the 
different models are fitted with an eight order polynomial 
in log t and log Z. The rms of the fitting is smaller than a 
few percent in all the bands. The polynomial fits provide 
a simple, analytic tool to calculate the derivatives ct,x and 
cz,x for an SSP model with given age and metallicity. In 
practise, we solved Eqs. [7J by evaluating the derivatives at 
an age of WGyr and solar metallicity. We verified that the 
results are very insensitive to the choice of such age and 
metallicity values. Moreover, using different stellar popula- 
tion models, as described in paper II, does not change at 
all the results. Eqs. [7] were solved in a least-squares sense, 
using the estimated values of ct,x and cz,x, and the values 
of b 



from Tab. [2] In order to estimate the uncertain- 
ties on the solutions, 6 log£/<51og Ejv and <51ogZ/<51ogEjv, 
we repeat the procedure by shifting the 6 Sjv x 's according 
to their uncertainties. 
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Figure 13. The peak value of the logi? e — alogcro distribution, 
computed by the bi-weight statistics, is plotted as a function of 
local galaxy density for the samples of ETGs in groups (filled 
circles) and the field sample (empty circle). The field sample is 
marked by the vertical dotted line. The best-fitting line to the 
data is plotted as a dashed line. The dot-dashed lines mark the 
±1 <r confidence levels on the best-fit relation. 




E N (gals/Mpc ) 

Figure 14. Same as Fig. 1131 but for the quantity logi? e ,r 
2 log (To, which is a proxy for galaxy dynamical mass. 



Fig. [15] plots the 6 Ejv x as a function of the effective 
wavelength of the filter X, from g through K. The solid 
curve is obtained by inserting the best-fitting solutions, 
8 log t/5 log E jv and S log Z/S log E jv , into Eq. [7] and com- 
puting the expected values of fc Ejv , x • We see that a combined 
variation of age and metallicity is able to match very well 
the value of 6„ „ in all wavebands. The dashed and dotted 
curves in the Figure are obtained with the same procedure 
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Figure 15. Variation of the FP intercept per decade in local 
galaxy density as a function of the effective wavelength of each 
filter (out of grizYJHK). The solid, dashed, and dotted curves 
are age+metallicity, pure age, and pure metallicity stellar pop- 
ulation models (see the text). Each waveband is indicated by a 
different colour, as shown by the labels beneath the error bars. 



by setting either 5 log Z/S log Ejv = or S log t/8 log Ejv = 
in Eq.0 They correspond to pure age and metallicity models 
of the variation of "c" with local density. The pure metal- 
licity model is clearly inconsistent with the observations. 
In fact, a difference in metallicity would tend to produce 
a negligible difference of mass-to-light ratios in K-band, as 
the NIR light is very insensitive to differences in metallic- 
ity, through line-blanketing. On the other hand, a pure age 
model is still able to match the observations, but with a 
large discrepancy, of ~ 2 a, in K-band. 

Fig. [16] plots the values of 8 log t/8 log Ejv vs. 
8 log Z/S log Ejv obtained from the procedure described 
above. The scatter in the plot (black dots) reflects the mea- 
surement uncertainties on the quantities b^ N x . As expected 
from the fact that 8 log Z/S log Ejv and 8 log t/8 log Ejv are 
estimated from the same system of equations (Eqs. [7]), 
the measurement errors on b s x imply a correlated vari- 
ation of 8 log Z/S log Ejv and <51ogt/<51ogEjv. A linear fit- 
ting of the 5 log Z/S log Ejv and 5 log t/8 log Ejv values 
(dashed line) gives a slope of —4.36 ± 0.05. Taking the 
mean and standard deviation of S log Z/S log Ejv 's and 
<51ogt/(5logEjv's, we obtain l S log Z/S log Ejv = -0.043 ± 
0.023 and (51ogi/51ogEjv = 0.048 ± 0.006, i.e. ETGs in 
higher density regions are older, by ~ 11% per decade in 
local density, and less metal-rich than galaxies at low den- 
sity. We notice that the value of (5 log Z/8\og Ejv is con- 
sistent with zero at about 2 a, implying that, as noticed 
above, a pure age model might be able by itself to ex- 
plain the observations. For a pure age model, the average 
value of 8 log t/8 log Ejv reduces to 0.04 ±0.004, rather than 
0.048 ± 0.006. 

In Sec. [7] we have reported some evidence that the vari- 
ation of "c" per decade in local density is stronger for galax- 
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Figure 16. Variations of age and metallicity per decade in lo- 
cal galaxy density, 8 log t/S log Sjv and 8 log Z/8 log Sjv, implied 
by the variation of the FP offset with log Ejv in the different 
wavebands. The best solutions of Eq. [7] are marked by the red 
circle. The black points are the values of 8 log t/8 log Ejv and 
8 log Z/ 8 log Sjv obtained by shifting the x according to the 
corresponding uncertainties. The dashed line is obtained by an 
ordinary least-squares fitting procedure, taking the arithmetic av- 
erage of the results obtained by changing the role of the dependent 
variable in the fitting. 

ies residing in more massive (relative to poor) clusters. As 
shown in Fig. 1111 for M group >2 x 10 14 Mq , the absolute value 
of 6 Sjv is about twice larger than that of —0.0097 reported 
in Tab. [2] For a pure age model, from Eq. [7] one would in- 
fer that, for more massive clusters, the variation in age per 
decade of local galaxy density is around twice larger than 
that reported above, i.e. Slog t/8 log Ejv ~ 0.08. 



9 SLOPES OF THE FP AS A FUNCTION OF 
THE ENVIRONMENT 

So far we have analysed the offset of the FP, interpreting 
it as a difference of the average mass-to-light ratio of ETGs 
residing in different environments. In this section, we study 
the dependence of "a" and "6" themselves on environment. 
In Sec. 19.11 we analyse the slopes "a" and "b" in the r band, 
taking advantage of the large size of the optical ETG's sam- 
ple. In Sec. 19.21 we derive the waveband dependence of FP 
slopes. 



9.1 The slopes of the FP in r band 

We start comparing the slopes of the FP for the r-band 
samples of field and group galaxies (Sec. I5.2"jl . To this aim, 
we have to account for two different effects, (i) field galax- 
ies have on average lower mass-to-light ratios than those in 
groups (see Sec. [8}, and (ii) ETGs are characterized by a 
luminosity segregation, for which brighter ( i.e. higher L) 
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galaxies mostly inhabit the denser cluster regions. More- 
over, not only L, but also galaxy radii are expected to change 
with environment, as they correlate with position within 
cluster (e.g. Cypriano et al. 2006). Point (i) implies that, at 
a given mass, an ETG in a lower density region is on average 
brighter than a galaxy with the same mass, but in a denser 
environment. Moreover, even removing this effect, the rel- 
ative distribution of galaxy luminosities and mean surface 
brightnesses might depend on environment (point ii) . To ac- 
count for points (i) and (ii), we proceed as follows. First, we 
fix the slopes of the FP, as in Sec. 17.11 and derive the inter- 
cepts of the FP for the field and group samples. From Eq. SI 
the average luminosity difference between the two samples 
is given by SE = -(8c) /b = -0.077 ± 0.007 mag. Then, we 
shift the < n > values of galaxies in the field sample by 
—SE, and bin both the field and group samples with respect 
to total magnitude and mean surface brightness. For each 
bin, we randomly extract the same number of field and group 
galaxies. This procedure, described in detail in paper II (see 
the app. B), allows us to extract two subsamples of field and 
group galaxies having the same distribution in the space of 
the effective parameters. By construction, both subsamples 
consist of the same number of galaxies, i.e. N = 11, 024 out 
of 11, 824 and 16, 717 ETGs in the field and groups, respec- 
tively. It follows that any difference in the FP slopes for these 
two subsamples arises from the different ways velocity dis- 
persions are related to the effective parameters, rather than 
from any spurious ( "geometric") effect due to different rela- 
tive fraction of galaxies in different regions of the parameter 
space. In other words, we are "removing" the dependence 
of galaxy luminosities and radii on environment, in order to 
measuring the genuine dependence of the FP on environ- 
ment. 

Fig. [T7] plots the "a" and "b" values for the field and 
group subsamples. We find a significant difference between 
the FP of field and group galaxies, with field galaxies having 
higher "a" and lower "b" value. For the "a", the difference is 
significant at ~ 3 a. Notice that had we not homogenised the 
luminosity and < fi > c distributions of galaxies in the two 
samples, we would have found an almost negligible differ- 
ence in "a" (at ~ la), and a stronger difference in "6" (at 
~ 5 a). The dependence of galaxy luminosities and radii 
on environment plays an important role. If not accounted 
for, it would lead to underestimating the dependence of the 



coefficient 



on environment. To further test this result, 



we have implemented a procedure to wash out the depen- 
dence of FP on environment and measure only the effect 
that the environmental dependence of galaxy luminosities 
and radii have on FP coefficients. For a given galaxy pa- 
rameter (e.g. luminosity), we bin the sample of field+group 
galaxies with respect to that parameter, each bin including 
the same number of N — 100 galaxies. Then, we scramble 
the galaxy environments, preserving the number of objects 
residing in each environmental bin (i.e. field and groups). 
The resulting shuffled sample is split again into field and 
group galaxies. Since both galaxy luminosities and radii are 
known to be affected from environments (see above), we per- 
form the scrambling with respect to both L and r e . Also, we 
consider the case of scrambling < /i > c (as it reflects both 
L and r e ) and the remaining FP variable, logoo. The FP 
slopes of the scrambled samples of field and group galaxies 
are reported in Tab. [3] The coefficient "b" is fully consistent 
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Figure 17. Slopes "a" and "6" of the FP for field (black circle) 
and group (red cross) galaxies. The solid, dashed, and dot-dashed 
ellipses mark the 1, 2 and 3 a confidence contours on "a" and 
"6", respectively. Notice that both the field and group samples 
have the same distribution in the space of the effective parameters 
(after the average difference in mass-to-light ratio is accounted 
for), and include the same number of 11, 024 ETGs (see the text). 



between the scrambled samples of field and group galaxies. 
For the "a", we see that scrambling < fi > e and logoo, 
does not impact at all its dependence on environment. On 
the other hand, scrambling L and r e has the net effect of 
decreasing (by ~ 5%) the "a" for field (relative to group) 
galaxies. This result is fully consistent with our conclusion 
above. Field galaxies have higher "a" than group galaxies. 
However, the dependence of L and r e on environment does 
not make this difference immediately detectable, as it tends 
to decrease "a" for the field sample. 

We derive the FP slopes of ETGs residing in groups 
with and without substructures. The presence of substruc- 
tures is characterized by a 3D A test (see Sec. 14. 2\ . About 
44% of groups show evidence of substructures. ETGs in the 
optical sample are first split in two subsamples according to 
the presence (N = 9133) or not (N = 7584) of substruc- 
tures in the parent groups. Notice that although less than 
half of the groups have substructures, there is a large num- 
ber of ETGs residing in groups with substructures. Then, 
proceeding as described above, we extract two subsamples 
of galaxies with the same distribution in the space of effec- 
tive parameters. These subsamples include 7, 716 galaxies 
each. No correction for the average mass-to-light ratio is ap- 
plied, as we obtain similar values of the FP offset for the 
two samples, i.e. c = -8.876 ±0.002 and c = -8.882 ±0.002 
for the case with and without substructures, respectively. 
Fig. 1 181 compares the FP slopes for the two "homogenised" 
subsamples. The corresponding values of "a" and "6" are 
fully consistent. 

To compare the FP slopes of galaxies residing in differ- 
ent bins of local density, we apply the same binning proce- 
dure used to characterise the trend of "c" with Ejv (Fig. [6j) . 
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Table 3. FP slopes of ETGs in groups and field obtained by scrambling galaxy 
environment at fixed size, surface brightness, velocity dispersion, and logarithmic 
luminosity. 



scrambling field groups 

a b a b 





1.350 ±0.019 


0.315 ± 0.001 


1.406 ±0.014 


0.316 ±0.001 


<M>c 


1.399 ±0.016 


0.315 ±0.001 


1.388 ±0.016 


0.316 ± 0.001 


log (T 


1.401 ± 0.017 


0.315 ±0.001 


1.391 ±0.014 


0.316 ± 0.001 


logL 


1.353 ±0.020 


0.314 ±0.001 


1.402 ±0.014 


0.316 ± 0.001 
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Figure 18. Slopes "a" and "5" of the FP for galaxies residing in 
groups with (red) and without (black) substructures. The solid 
and dashed ellipses mark the 1 and 2 <r confidence contours on 
"a" and "b", respectively. 



To derive the slopes of the FP in each bin, we proceed in a 
similar way as in the comparison of field and group galax- 
ies, i.e. we account for (i) the different average mass-to-light 
ratio of ETGs in different environmental bins, and (ii) dif- 
ferences in the relative distribution of galaxies in the space 
of effective parameters. First, we estimate the average lumi- 
nosity difference of galaxies residing in different bins with 
respect to the field sample, using the linear fit of "c" ver- 
sus logEjv (Eq. [5j. For a given density bin, i, the average 
luminosity difference is given by 



5Ei = — 



X (logEjv - logE/ieid) 



(8) 



where ~Ef ie id is the local density value assigned to field galax- 
ies (Sec. 17. ip . We then shift the <n> e values of the entire 
optical sample of ETGs by 8Ei , and extract a subsample of 
galaxies with the same luminosity and </j,> c distributions 
as galaxies in the i-th Ejv bin. We compute the relative vari- 
ations, 8a/ a and Sb/b, of the FP slopes of the entire sample 
after these selections are applied. The factors 6a/a and Sb/b 
are applied to the FP slopes of galaxies in the given local 
density bin. By construction, the correction vanishes for the 
field sample, while it become more and more important as 



Ejv increases. Fig. [19] plots the r-band slopes of the FP as 
a function of local galaxy density. In order to characterise 
the trend of "a" and "b" with Ejv, linear fits of both slopes 
with respect to logEjv are performed. Notice that the field 
sample is not included in the fit. 

In general, one can notice that the FP slopes exhibit 
only a weak dependence on environment. Before applying 
any correction, the "a" does not show any variation with 
Ejv, with the slope of the linear regression (0.00 ± 0.03; see 
Fig. I19|l being fully consistent with zero. After the correc- 
tions are applied, we are able to detect a 2.5 a tendency 
of "a" to decrease from low to high density environments. 
One can notice that, in agreement with this trend, the cor- 
rected value of "a" for the field sample is larger than all 
the slope values obtained for the higher density subsam- 
ples (Ejv > Igals/Mpc 2 ). This is consistent with what 
found above, when comparing the FP slopes of field and 
group galaxies (Fig. I17|l . Notice that the values of "a" and 
"b" in Fig. [T7] refer to a different sample of field galax- 
ies with respect to that considered in Fig. 1191 In fact, for 
Fig. 1171 two subsamples of field and group galaxies are ex- 
tracted, to match the corresponding distributions of galaxies 
in the space of structural parameters. For the comparison in 
Fig. 1191 the FP slopes are derived by using all galaxies in the 
field sample. This explains the difference (in absolute value) 
of the values of "a" and "6" for field galaxies between Fig. [17] 
and Fig. 1191 Despite that, the relative difference of FP co- 
efficients between the field and group samples are detected 
in both cases. The coefficient "6" tends to increase with lo- 
cal galaxy density. The correction procedure shows that the 
tendency is significant at ~ 3a, and is mostly due to galaxies 
in higher density regions (Ejv>10 gals Mpc~ 2 ). As shown in 
Fig. 1201 the trends of "a" and "6" with environment be- 
come weaker when plotting their value as a function of the 
normalised cluster-centric distance R/R200- In this case, the 
corrections are estimated with the same approach described 
above, using the fictitious value of the cluster-centric dis- 
tance of the field sample (Rfieid/R20o) and the linear fit of 
"c" versus R/R200 (see Fig. [7J. We find that the "6" in- 
creases from the centre to the periphery of galaxy groups, 
as expected from the variation of "b" with local galaxy den- 
sity. For the "a", we do not detect any significant variation 
with R/Rwo- 

Fig. [^Uplots the FP slopes in r band as a function of the 
mass of the parent galaxy groups. As for Ejv and R/R200, 
we plot the values of "a" and "b" with and without the 
corrections for the effects (i) and (ii) mentioned above. We 
find that the coefficient "a" tends to decrease (at ~ 2 a) 
for galaxies residing in more massive clusters, while the "6" 
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Figure 19. The slopes "a" and "&" of the FP are plotted as 
a function of local galaxy density, Sjv- Error bars denote 1 a 
uncertainties. Red colour shows the case where the slopes are 
corrected for the environmental dependence of the average mass- 
to-light ratio of ETGs and luminosity segregation of galaxies in 
different environments (see the text). The dashed lines are the 
linear least-squares fits to the data points. The squares in the left 
part of the plots are the slope's values for the field sample, as also 
indicated by the vertical dotted line and the arrow. 



shows an opposite and more significant (at ~ 3 a) trend, 
increasing for higher M group . A similar result is obtained 
when characterising the global environment in terms of the 
"global" projected mass density, M group / -R|oo) of the parent 
galaxy groups. The "a" ("&") exhibits only a weak trend to 
decrease (increase) with M group / R\ 00 (Fig. [22}. As shown by 
the slope values of the linear fits, the trends of "a" and "b" 
with "global" density are only marginally significant, at less 
than 2 a. On the contrary, the dependence on local density 
seems to be more significant, as shown by the uncertainties 
on the slope values reported in Fig. 1191 This implies that the 
local rather than the global environment is the main driver 
of the environmental dependence of the slopes of the FP. 
The same result has been shown to hold for the offset of the 
FP (Sec. ED- 

As for the "c", we also analyse how the slopes of the 
FP depend on local galaxy density in different bins of mass 
of the parent galaxy groups. Figs. [23] and [24] plot the "a" 
and "b" as a function of log Ejv, in different bins of M group . 
The binning is performed in the same way as for Fig. 1101 
The slope values are corrected as in Fig. 1191 Fig. [23] does 
not show any significant dependence of "a" on local galaxy 
density for all bins of M group . This result is not necessar- 
ily inconsistent with the trend of "a" with Ejv shown in 
Fig. [T5] (upper panel) , as binning the sample with respect 
to both Mgroup and Ejv decreases the sample size in each 
bin, possibly washing out the (2.5 a) environmental trend in 
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Figure 20. The same as Fig.[l9]but plotting the "a" and "b" as 

a function of the normalised distance to the centre of the parent 
groups where galaxies reside. 
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Figure 21. The same as Fig. ll9l but plotting the "a" and "6" as 
a function of mass of parent groups where galaxies reside. 
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a function of "global" surface mass density, M group / R\ 00 , °f * ne 
parent groups where galaxies reside. 



Fig. 1 191 Fig. 1241 shows that "6" does not depend significantly 
on group mass for low mass groups, while it increases signif- 
icantly with log E at in the high mass regime. In the bin of 
highest Mgroup (5.8 x 1O 14 M ), the slope of the "6"-log£jv 
linear fit is significantly positive, at more than 4 a, for the 
case where no correction is applied, and at more than 2.5 a, 
when the correction is performed. We conclude that the de- 
pendence of "6" on local density (Fig. I19|) is mainly driven 
by galaxies residing in most massive clusters. 



9.2 The slopes of the FP from g through K 

The environmental analysis of the waveband dependence of 
the FP involves essentially two aspects. We derive the rela- 
tive variation of FP slopes, from g through K, and compare 
such variation among different environmental bins. Since 
the waveband variation of FP slopes informs on how stel- 
lar population properties (i.e. age and metallicity) vary with 
galaxy mass (see paper II for details), this allows us to anal- 
yse how such stellar population variation is affected by the 
environment where galaxies reside. Moreover, as the NIR 
light closely follows the stellar mass distribution in galax- 
ies and the contribution of stellar populations to the tilt 
of the FP in K-band vanishes (see paper II), the NIR (K- 
band) FP informs on how the ratio of dynamical to stellar 
mass (and/or the non-homology) of ETGs vary with mass, 
and how this variation is related to the environment. The 
results presented below are discussed, with regard to both 
aspects, in Sec. 10.2. 

In order to derive the FP coefficients in different envi- 
ronments, we use the optical+NIR sample of ETGs, binning 
galaxies with respect to different environmental properties, 



and deriving the FP slopes in each bin. In all cases, the bin- 
ning is performed as for the environmental analysis of the 
FP intercept from g through K (Sec. I7.2[l . We correct the 
FP slopes for the same effects mentioned in Sec. 19.11 i.e. 
environmental differences in the average mass-to-light and 
luminosity segregation. To this effect, we follow the same 
recipes as for the analysis of FP slopes in r-band (Sec. 19.1]) . 

Fig. [25] plots the slopes of the FP from g through K 
for galaxies residing in groups and in the field. The opti- 
cal+NIR samples of group and field galaxies consist of 2, 185 
and 1, 359 galaxies, respectively. The upper panel of the Fig- 
ure shows the FP slopes for these samples: we only detect 
a significant difference in the value of b between field and 
group galaxies, with the latter having lower "6" in all wave- 
bands. However, to perform a proper comparison, we have to 
account for differences in mass-to-light ratio and luminosity 
segregation. To this effect, we proceed as for the comparison 
of field and group galaxies in the r band (see Sec. 19. ip . by 
extracting two subsamples of galaxies with the same distri- 
butions in the space of structural parameters. This selection 
leads to field and group sub-samples with 1, 264 galaxies 
each. The corresponding FP slopes, from g through K axe 
shown in the lower panel of Fig. 1251 We are now able to 
detect also a difference in the coefficient "a". The "a" turns 
out to be smaller for group, relative to field, galaxies in all 
wavebands, in agreement with the result for the r-band sam- 
ples (Tig. ll7|) . In order to measure the variation of FP slopes 
from g through K, we assume a linear relation between the 
"a" ("6") and the effective wavelength (A) of each waveband, 
i.e. 



a(A) = p a + (?alogA 
6(A) = p b + q b log A, 



(9) 
(10) 



where p a {pb) and q a (qb) are the intercept and slope of 
the linear relation. The values of p a and q a {pb and qb) are 
derived by an ordinary least-squares fit of "a" ("6") ver- 
sus log A. The uncertainties on p a and q a {pb and qb) are 
estimated from those on "a" and "6", accounting for the 
correlated errors of FP slopes among different wavebands Q. 
We find that, for all the subsamples of ETGs analysed in the 
present paper, Eqs. l9l and fTDI are able to describe the trends 
of "a" and "b" with waveband with an accuracy of a few 
percent. They provide simple, empirical relations to quan- 
tify the waveband variation of FP slopes. For both field and 
group ETGs (lower panel of Fig. [25}, we use Eqs. [9l and[T0l 
to calculate the relative variation of "a" and "6" from g 
through K, i.e. 



da 



g->K 



fib. 



g->K = 



a(\ K ) - a{Xg) 
a(X g ) 

b(X K ) - b(X g ) 



(ii) 

(12) 



4 Since we analyse the same sample of ETGs and velocity disper- 
sions are the same in all different wavebands, the uncertainties on 
"a" and "6" from g through K are correlated. To estimate the 
errors on p a and q a (p b and q b ), we randomly shift the values of 
"a" ("6") from g through K, and repeat the fitting. The shifts are 
computed by taking into account the correlation of uncertainties. 
Had we not accounted for it, the errors on p a and q a (pb and q b ) 
would have been significantly overestimated, by almost a factor 
of two. 
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Figure 23. Dependence of the slope "a" of the FP as a function of local density, Sjv, in different bins of mass of the parent galaxy 
group. Red colour refers to the values of the slope after correcting for differences in the FP intercept and distribution in the space of 
structural parameters (see the text). Error bars denote 1 a uncertainties. For each panel, the corresponding median value of the mass of 
the parent galaxy group, < M gr0 up >, is reported in the upper-right. The < M gr0 up > increases from left to right, and top to bottom. 
The values of "a" for the field sample are reported in each panel, and are marked by the vertical dashed line. The dashed lines are best-fit 
linear relations to the data points, with the corresponding slope values being reported in the lower-left part of each panel. Notice that 
the S jv bins are the same as in Fig. 1101 



where A 9 and Xk are the effective wavelengths of the g 
and K passbands, respectively. The values of Sa g - > K and 
Sbg-yK define a segment in the "a"-"b" plane, whose size 
quantifies the amount of variation of FP slopes from the op- 
tical to NIR. The segments corresponding to the subsamples 
of field and group ETGs, as well as the values of 5a 9 _>K 
and 8b g - > K, are shown in the lower panel of Fig. 1251 The un- 
certainty on Scig-yK (Sbg-ypc) is obtained from those on p a 
and q a (pi, and q ). In general, the variations of "a" and "6" 
from g through K are similar for both samples. For the "b", 
field galaxies tend to have a marginally steeper variation, 
with Sb g - >K = 2.5±0.7%, compared to 8b g - >K = 1.1±0.7% 
for group galaxies. More in general, one can notice that while 
the "a" is always smaller for group relative to field galax- 
ies, the difference in "b" between the two samples tends to 
vanish when moving from "g" to "K" (Sec. 19. 1|) . To fur- 
ther quantify these differences, we take the average of the 



g- and r-band coefficients, and compare them to the mean 
values of "a" and "b" in the JHK (NIR) wavebands. For 
the field (group) sample, the "b" changes from 0.302 ±0.003 
(0.312±0.002) in the optical to 0.311 ±0.003 (0.316±0.003) 
in the NIR. Thus, the difference in "b" between the field 
and group samples is significantly in the optical, while be- 
comes consistent with zero in the NIR. For the field (group) 
sample, the "a" changes from 1.35 ± 0.06 (1.23 ± 0.04) in 
the optical to 1.53 ± 0.06 (1.42 ± 0.05) in the NIR. Thus, 
the difference in "a" between the field and group samples 
is essentially independent of waveband, with the "a" being 
smaller (at ~ 2 a), for the group relative to field sample in 
both the optical and NIR. 

Fig.[26]plots the FP slopes from g through K for galax- 
ies residing in groups with and without substructures. We 
do not find any significant difference, for all wavebands, be- 
tween the two samples. 



Environmental Dependence of the FP 21 



d 



Q 



d 



i*5 



d 



tttt] — i — i i nun — 

field <m„, 



.""l I 

field <i»„, 



ii| — t-ttt 

, >-0.2x10 ,4 M, 



TTJTTj 1 — TTTTTTT] 

, ifield <M_, 



i 1 1 1 ni| — i i n 
>=0.3x10 14 M, 



-I 



ih ii| i iii n ii| i i n ii n l i i i_mhh | i i m iii i i| i ii _mh ii| i :iii n ii| i i M i ni ' ' U 



,„>=0,1«10' , M, 



slope; 



iC*- — 



).006±0,005 



slope=«.001 ±0.004 



slope— .002 ±0 . 005 
slope=-0.001 ±0.005 



slope«0.005±0l005 
■ alopB=0.001 ±0.005 



ifield <m^ u ^>-o.5xio 14 m, 



'field <M, rmip ^-o.7xio r ,4 M, 



WW 



slope -0.002±0.005 

slope--0.000±0.005 



ifield <m ^-i.ixio'V 



• ir. 



slope- 0.006*0.005 
slope=- 0.00 1 ±0.006 



h*h| — I i n i m | — i i i m il ] — H-yfHoj — i i — i i m i n i] — i-i-h* 



slope«0.006±0.005 
slope- -0.001 ±0.005 
_luJ 1 — I I nu l l 1 — I 1 — i_l 



field <m >=i.9xio H M ( , Ifield <w ansii(! >=5.Bxio u w, 0.1 1 10 



S — 



/\ slope-0.009±0.005 

slope-0.003±0.005 



3 Op 



I I I N (gals/Mpc 2 ) 
A 



pe-0fO19±O.OO5 



3lope-sGU)13±0.005 



""I I ' ■ l i mi | I r lihI | i i t fliij i / i i Tijl | ill 



0.1 



1 



10 



0.1 



1 



10 



Z N (gals/Mpc 2 ) I N (gals/Mpc 2 ) 



Figure 24. The same as Fig. [23] but for the coefficient "6" of the FP. 



Figs.[27land[28lplot the slopes of the FP from g through 
K in different bins of local galaxy density, and mass of the 
parent galaxy groups. In each environmental bin, we apply 
suitable correction factors to the slope values, in order to ac- 
count for environmental differences in the average mass-to- 
light and luminosity segregation. The corrections are com- 
puted as in the analysis of FP slopes in r-band (Sec. l9.ll) . For 
all wavebands, we apply the same corrections. This is justi- 
fied by the fact that (a) we always analyse the same sample 
of ETGs in the different wavebands, and (b) as shown in 
Sec. I7.2I the amount of environmental variation of the FP 
offset is very similar in all wavebands, implying similar dif- 
ferences in the corresponding average mass-to-light ratios. 
We find that the waveband dependence of the FP is very 
similar in all wavebands. The "a" tends to increase with 
wavelength in all environments. We measure the amount of 
end-to-end variation of "a" and "b", across the available 
wavelength baseline, by computing the quantities 5a g _>K 
and <5& 3 ->/f(see above). The value of 8a g ->K spans the 
range ~ 11% to ~ 23%, with a weighted-mean value of 16%. 
All values of 8a g - > K are consistent within the errors with 
that of 16 ± 3% measured for the field sample. The "6" does 



not change significantly with waveband, with a marginal ten- 
dency to increase (Sb g ->K> 0) by 1-3% from g through K 
in all the environments. Notice that the quantity Sbg-yK 
is positive in all bins of Ejv and M group with the excep- 
tion of the the bin with highest value of parent group mass 
{Mgroup = 3.05 x 10 14 M Q , see Fig. [28]). This is consistent 
with what shown in Fig. 1251 i.e. the "6" increases by a few 
percent with waveband for the field sample, while it is con- 
stant for galaxies in groups. Fig. [28] would also imply that 
this different behaviour is due to galaxies residing in massive 
clusters. The plots of "a" and "6" as a function of other en- 
vironmental quantities, such as the cluster-centric distance 
and the "global" mass surface density of the parent galaxy 
groups, do not add further information to the conclusions 
above, and are not shown here for brevity reasons. 

To summarise, we find that the tilt of the FP, as mea- 
sured by the coefficient "a", is different between field and 
group galaxies, with group galaxies having lower "a". The 
difference is independent of the passband where the galaxy 
structural parameters are measured. For the "6", its varia- 
tion with waveband is smaller than ~ 3%. However, in this 
case, an important difference between field and group galax- 
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Figure 25. Slopes "a" and "ft" of the FP from g through K for galaxies residing in groups (empty circles) and in the field (filled circles). 
Different colours correspond to different wavebands, as shown in the lower-left corner of the upper plot. Ellipses denote 1 a confidence 
contours on "a" and "ft". The upper and lower panels correspond to the cases where the FP coefficients are either corrected or not for 
the difference in the average mass-to-light ratio of galaxies in the different environments and the luminosity segregation effect. The black 
arrows in the lower panel plot the amount of relative variation from g through K of the FP slopes for the two subsamples, as obtained 
by assuming a power-law relation between the "a" ("6") and the effective wavelength of a given passband (see the text). The relative 
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Figure 26. The same as the lower panel of Fig. l25l but for galaxies 
residing in groups with (filled circles) and without (empty circles) 
substructures. 



ies emerges, i.e. the "6" increases with waveband (by 2.5%, 
see the value of &b g - > K in Fig. I25[) . while for group galax- 
ies the variation with waveband is insignificant. This seems 
to be especially true for galaxies residing in most massive 
clusters (Fig. [28}, for which 56 9 _>k is negative. The impli- 
cations of these findings are discussed in Sec. 1101 



10 DISCUSSION 

In the present work, we have analysed the environmental 
dependence of the FP coefficients on the waveband at which 
structural parameters are measured. On what follows, we 
discuss the environmental dependence of the FP in the op- 
tical (Sec. 10.1) and using the combined optical+NIR (Sec. 
10.2). In Sec. 10.3 we examine the stellar content of ETGs 
in different environments. The most remarkable result we 
find here is the variation of the FP offset with local density. 
As shown in Sec.[S] this trend sets strong constraints on the 
variation of stellar population properties of ETGs with the 
environment. This point is discussed in Sec. 110.31 



10.1 Environmental dependence of the optical FP 
relation 

The first attempt to measure systematic differences in 
the FP relati on among different environment s has been 
performed by Ide Carvalho fc Diorgovsk (|l992T ) (hereafter 
dCD92; see also iDiorgovskil Ide Carvalho ) . who derived FP 
coefficients in the optical for a sample of 31 cluster (mostly 
Coma and Virgo) and 34 "field" galaxies in the low-redshift- 
Universe. Field galaxies included also objects residing in 
poor and loose groups. The authors found that the offset 
"c" of the FP for field galaxies is significantly different from 
that of galaxies in clusters. They interpreted qualitatively 
this difference as field galaxies being younger than those in 
clusters, consistent to what we find here. However, a direct 
and quantitative comparison of the values of "c" between 
our study and dCD92 is not possible because of the differ- 
ent definition of environment and the different techniques 
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Figure 27. Slopes "a" and "6" of the FP in different bins of local galaxy density. The average value of log Sjv in a given bin is reported 
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Figure 28. The same as Fig. [27] but for different bins of mass of parent galaxy groups. 
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used to estimate the FP intercept. Due to the uncertain- 
ties on FP coefficients, dCD92 did not detect any significant 
variation of FP slopes with environment. 

The first systematic study of the FP relation 
in clusters, at optical waveban ds, wa s performed 
by Ijoreensen. Franx. fc Kiaergaardl (|l996l ) (hereafter 
JFK96). They derived the FP coefficients for ETGs, in 
Gunn-r, residing in 11 nearby clusters (220 galaxies). They 
found no evidence for variation of FP coefficients with 
global cluster properties like richness, velocity dispersion 
and temperature of the intra-cluster gas. Also, they did not 
detect any correlation between residuals around the FP and 
the local environment, characterized as the cluster-centric 
dista nce and the projected cluster surfa c e den sity (see 
also lLa Barbera. Busarello. fc Capacciolil I2000T ). These 
results seem to contrast to what we find here, i.e. that 
both "c" and (to minor extent) "a" and "6" depend on the 
environment (Sees. [7] and [9l respectively). However, one 
has to consider that the number of galaxies per cluster in 
JFK96 is much smaller (larger uncertainties) than that in 
each local density bin analysed here. For instance, looking 
at tab. 3 of JFK96, one can see that their typical error 
on "b" varies from ~ 0.02 to ~ 0.06, among the different 
clusters, while in our case this error amounts on average to 
~0.006. The environmental variation of "6" we find ranges 
from ~0.02 to ~ 0.04 (see the panel representing the more 
massive systems in Fig. [24}, comparable to the quoted 
uncertainties on "b" in JFK96. This explains why we are 
able to detect such a small difference while JFK96 cannot. 

iBernardi et al.l |2006l ) (hereafter B06) analysed the off- 
set of the FP relation of ETGs in low and high density re- 
gions (bins), using a large sample of ETGs from SDSS. The 
authors defined the environme nt using the C4 clu ster cata- 
logue (based on SDSS-DR2; see lMiller et~ai1l2005h . Galaxies 
were assigned to the high density bin using a cut in angular 
separation and radial velocity around each cluster, rather 
than assigning membership to each object as we do here 
(Sec. [5}. Galaxies in the low density bin were defined us- 
ing a cutoff in the maximum allowed distance for a " field" 
galaxy to the n-th companion galaxy. This criterion was in- 
troduced by B06 in order to remove the contamination in 
the field sample from galaxies residing in poor and loose 
groups, which are likely not included in the C4 catalogue. 
On the contrary, the FoF catalogue we use here also includes 
poor galaxy groups, and we define the field sample accord- 
ing to the 3D distance of galaxies to the nearest group. De- 
spite of these differences and the different sources of struc- 
tural parameters used to derive the FP in the present study 
(i.e. 2DPHOT, see paper I) and B06 (i.e. the SDSS Photo 
pipeline), the difference of the FP intercept among field and 
group galaxies between our study and that of B06 is very 
consistent. In fact, the latter find that the average difference 
of < jj, > c in r band, as derived from the FP, amounts to 
0.075 ± 0.008, with galaxies in denser regions having fainter 
surface brightness than their counterparts in the field. On 
the other hand, we find -Sc/b = 0.077 ± 0.007 mag (see 
Sec. 19. 1|) . which is fully consistent with the B06 result. We 
address the implications of this result in Sec. 10.3. 

A systematic study of the environmental effects on 
F P coefficients has been recently performed, in the optical, 
bv lD'Onofrio et all d2008h (hereafter D08), in the framework 
of the WINGS project (|Fasano et aljfeood ). The sample of 



D08 consists of 1, 550 early- type galaxies in 59 massive clus- 
ters, in the low-redshift-Universe (0.04 < z < 0.07). Hence, 
they analyse a similar redshift range to that considered by 
our study, but with a different range in mass of the parent 
systems where galaxies reside. D08 detected no significant 
variation of FP coefficients with global cluster properties, 
such as richness, -R200, and velocity dispersion. On the con- 
trary, we find evidence that the offset "c" and (to minor 
extent) the slope "6" of the FP change systematically with 
the global environment, with "6" ("c") increasing (decreas- 
ing) with respect to the parent group mass, M group (see 
Figs. [5] and 1210 . We notice that these trends are proba- 
bly detectable only because of the larger range in M grouv 
spanned by our group/cluster catalogue (see Fig. 3) in com- 
parison with theirs (400 < a < 1400 km/s, with < a > = 
800 km/s). 

Moreover, as discussed in Sec. [7] the trend with global 
environment can be entirely explained by the dependence of 
FP coefficients on local density. D08 found a strong variation 
of FP coefficients with the local environment. Both "a" and 
"6" were found to increase, and "c" was found to decrease 
with respect to the (projected) local galaxy density. The 
trends of "6" and "c" from D08 are qualitatively consistent 
with those we find here (see Secs. l7,ll and l9.1[l . However, after 
correcting for systematic effects (see red circles in Fig. I19[) , 
we detect a 2.5<r indication that the "a" decreases from the 
low to the higher density regions, in contrast with D08. Con- 
sidering ETGs in the most massive groups (see Fig. 123ft . we 
still do not find any strong variation with local environment, 
as in D08. A possible explanation for this discrepancy can 
be the fact that we have corrected the FP slopes of different 
samples of ETGs for biases due to the environmental de- 
pendence of the average mass-to-light ratio of galaxies and 
the distribution of ETGs in the space of effective parame- 
ters, while D08 did not account for either effects. We also 
notice that the variation of "6" with local environment from 
D08 is much larger than that detected here. From Fig. 1191 
we see that the end-to-end variation of "fc" with respect to 
local density amounts to less than 0.02, while D08 find a 
variation of ~ 0.04 (see their fig. 15). This discrepancy can 
be explained by the different mass regime covered by our 
cluster sample and that of D08, and by the fact that we 
find the variation of "6" with local environment to depend 
on the mass of parent galaxy groups. For groups as massive 
as ~ 5.8 x 10 14 Mq, we find an end-to-end variation of "6" 
of ~ 0.04 when none of the above corrections is applied (in 
agreement with D08). This variation reduces to ~ 0.02 after 
the corrections are performed. 

10.2 Environmental dependence of the FP from g 
through K 

For what concerns the waveband variation of FP coefficients, 
we find that the "a" increases, by ~ 16%, from g through 
K, independent of the environment, while the variation of 
"6" with passband is smaller than ~ 3%. For field galax- 
ies, (see Fig. 1250 . there is a 3 a evidence that the "6" in- 
creases, by ~ 2.5%, from the optical to NIR, with the dif- 
ference in "6" between field and group galaxies vanishing 
in the NIR. For group galaxies, in particular those resid- 
ing in most massive systems (Fig. I28[) , no variation of "6" 
with waveband is detected. As discussed in paper II, the 
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waveband dependence of FP slopes informs on the varia- 
tion of stellar population properties with mass. For the en- 
tire SPIDER sample, the relative variation of "a" and "6" 
from g through K implies that ETGs, as a whole, have es- 
sentially coeval stellar populations, with metallicity being 
larger in more massive galaxies. As shown from eq. 7 of 
paper II, a smaller value of "6" (as that we find here for 
field galaxies) implies a steeper variation of the mass-to- 
light ratio of ETGs along the galaxy (mass) sequence. As 
the "6" increases from the optical to NIR for field galaxies, 
the corresponding mass-to-light ratio variation is related to 
the way stellar population properties change with mass. In 
other terms, the environmental variation of "6" might im- 
ply that group galaxies exhibit a shallower relation between 
stellar population parameters (i.e. age and metallicity) and 
m ass. Interestingly, this is consistent with the recent findings 
of lPasauali et al'l l|201Ch (hereafter PGFIO), who found that 
galaxies in denser environments have shallower age- and 
metallicity-mass relations. This is also a natural expecta- 
tion of hierarchical models of galaxy formation, as galaxies in 
denser environments are those accreted in the group environ- 
ment at earlier times, having their star formation quenched 
earlier through environmental-driven effects (e.g. strangu- 
lation), implying a shallower relation between (luminosity- 
weighted) age and stellar mass. A flatter metallicity-mass 
relation in denser environment might also result from the 
effect of tidal st ripping on the more boun d stellar material 
in galaxies (e.g. iKlimentowski et al.l 120091 ) as discussed by 
PGFIO. 

On the other hand, the difference in "a" between field 
and group galaxies is the same at both optical and NIR 
wavebands. Group galaxies have a larger FP tilt, as mea- 
sured by the coefficient "a", likely because of a waveband- 
independent origin. A possible explanation is a different vari- 
ation of dark-matter fraction and/or non- homology along 
the mass sequence of ETGs between galaxies in different en- 
vironments. For instance, group galaxies might have a larger 
tilt because of a steeper relation between dark-matter frac- 
tion and mass, due to (i) tidal interactions between galaxy 
dark-matter halos and parent group halos, i.e. lower, rel- 
ative to higher, mass galaxies are stripped more of their 
dark-matter halo as they fall into the cluster potential well; 
(ii) central galaxies have higher dark-matter fractions, as 
they acquire part of the dark-matter halo of the in-falling 
satellites. The first scenario might result from the fact that 
the total time span, during which galaxies are stripped, is 
large r for low mass galaxies (see e.g. eq. 9 of lCattaneo et aD 
l2010h . However, tidal forces are stronger for more massive 
galaxies, as a s hort orbital d ecay time for L > L* galax- 
ies is expected |Barnedlf989h . The scenario (ii) might be 
more plausible, provided that stripped galaxies have mostly 
late-type morphology. This would not change significantly 
the peak value of dynamical mass of ETGs as a function 
of environment, consistent with Fig. 1141 and increase only 
the dark-matter fraction for higher mass systems, implying 
a smaller tilt of the FP relation. Notice that a steeper re- 
lation between stellar to dynamical mass ratio and mass in 
ETGs for group (rel ative to field) galaxie s has r ecently been 
suggested also from lShankar fc Bernardil (|2009h . as a result 
of the wider range of luminosity-weighted ages character- 
izing cluster relative to field galaxies. These points will be 
further analysed in a forthcoming contribution in this series, 



using estimates of dynamical and stellar mass for ETGs and 
analysing their dependence on environment. 

It is interesting to notice that the NIR FP, and its envi- 
ronmental dependence, can be connected to the s tellar mass 
FP recently derived by iHvde fc Bernardil (|2009h . The stel- 
lar mass FP is obtained by replacing luminosity (i.e. sur- 
face brightness) in the FP equation with stellar mass. Since 
the K-band light closely follows the stellar mass distribu- 
tion in galaxies and the contribution of stellar populations 
to the tilt of the FP vanishes in the NIR (see paper II), 
the K-band FP sh ould essentially co i ncide with the stellar 
mass FP. In fact, IHvde fc Bernardil (2009) report a value 
of a ~ 1.54 ± 0.02 for the stellar mass FP, while we mea- 
sure a ~ 1.55 ± 0.02 for the K-band FP. This reinforces 
the fact that the environmental difference of FP coefficients 
we detect in the NIR bands, among different environments 
(see above), implies an environmental difference in the re- 
lation between the ratio of dynamical-to-stellar mass and 
galaxy mass. Given the excellent agreement of the K-band 
and stellar mass FP's, one might actually use either rela- 
tions to derive ETG's stellar masses, and test how realistic 
are the estimates from stellar population models. 

10.3 Stellar content of ETGs in different 
environments 

The environmental variation of the FP intercept, '"c", im- 
plies that the average mass-to-light ratio of ETGs increases 
from low to high density regions. As seen in Eq. this vari- 
ation can be interpreted as a difference in the stellar content 
of ETGs, at fixed mass, provided that the ratio of dynam- 
ical to stellar mass, M/M*, does not increase, on average, 
with local density. There are several arguments against such 
an increase. First, Fig. [H] shows that the peak value of the 
quantity Oq x R e - a. proxy for dynamical mass - does not 
change with the environment. So, the M/M* can increase 
with local density, only if stellar mass, at fixed M, decreases 
with Eat. However, as galaxies are accreted into groups and 
clusters, they are expected to be stripped off of their dark- 
matter halos, resulting, eventually, into a larger M*, for fixed 
M, at high density. These arguments make plausible our as- 
sumption of constant M/M+, that allows us to use the offset 
of the FP as a tool to analyse the variation of stellar popula- 
tion properties as a function of en vironment. It is i nteres ting 
to note that, as shown in fig. 8 of iBernardi et all (|2006h . at 
fixed luminosity, L, dynamical mass does not depend on en- 
vironment, while the g — r colour (which is a good proxy 
for M*/L) does. This implies that, at fixed L, the M/M* 
should change with environment. We notice that this is not 
inconsistent with the above assumption of constant M/M*, 
as we are assuming that the M/M* is constant at fixed M 
(not L). 

We find a significant difference of the FP offset, in r 
band, between field and group galaxies, fully consistent with 
that measured from B06 (see previous section). Using spec- 
tral abundance indicators, B06 concluded that the difference 
in "c" is due to a difference in age between field and clus- 
ter ETGs, with field galaxies being younger by ~ 1 Gyr. 
Using the information encoded in the FP offset at different 
wavebands, we show that the variation of the FP intercept 
with local galaxy density implies a variation in logarith- 
mic (luminosity-weighted) age per decade of local density 
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of <51ogt/c51og Ejv = 0.048 ± 0.006. From Eq. \5\ taking into 
account the difference of "c" between field and group galax- 
ies (see previous section), the value of 8 log t/S log Ejv trans- 
lates to a relative difference in age of ~ 18 ± 2% between 
the field and group environments. For a galaxy formation 
epoch of 10 Gyr, this implies an age absolute difference of 
St ~ 1.8±0.2 Gyr, larger than what found by B06. The find- 
ing that field ETGs are younger than those in high density 
environments is qual itatively consis t ent w ith some previous 
studies (see Sec. [TJ. I Clemens et al.l (120091 ) found that field 
ETGs are ~ 2 Gyr younger than their cluster counterparts, 
in agreement with the value of St we measure here. The same 
amoun t of difference in age was also found bv lThomas et all 
(2005). Similar results, but base d on smaller samples o f 
ETGs, were also obtained b y fe.g.)|Kuntschner et al] ( 
Terlevich fc Fo rbes (2002 ); ISanchez-Blazquez et alj ( 
de La Rosa et alj (120071). Our r esults are also consistent 
with those of ICooper et all ()2010h . who found that, at given 
colour and luminosity (and hence stellar mass), galaxies with 
older stellar populations favour regions of higher overdensity. 
We notice that all these studies have been based on the anal- 
ysis of spectral features (i.e. line indices) of ETGs. Hence, 
the inferred differences in age refer to the inner galaxy re- 
gion, typically inside one effective radius. As already noticed 
in paper II, the information provided by the waveband de- 
pendence of the FP is more related to the global properties of 
the light distribution in galaxies (as measured by the struc- 
tural parameters). In this regard, it is not affected by ra- 
dial population gradient s in ETGs and the i r poss ible depen- 
dence on environm ent l)La Barbera et aD I2005T ) . Recently, 
iRogers et al. (2010) have performed a Principal Component 
Analysis of the SDSS spectra of ETGs residing in groups 
and clusters, spanning a wide range in mass (similar to that 
of our FoF catalogue). They find galaxies populating the 



lowest mass halos (M 9 



10 Mq) to be younger, by 



~ 1 Gyr, than those in most massive clusters. Using the lin- 
ear fit in Fig. [8] we estimate a variation of the FP offset, over 
the entire range of parent group mass, of ~ —0.015. From 
the measured values of fc Sjv and <51ogt/<5Ejv (see Sec. |8]), 
this difference translates to a difference in age of ~ 12%, 
i.e. ~ 1.2 Gyr (for a form ation epoch of 10 Gyr), consistent 
with lRoeers et alj l|20ld ). 

Recently, iThomas et aD (|20Qg| ) have analysed a large 
sample of morphologically selected ETGs from the SDSS, 
deriving stellar population properties, i.e. age, metallicity, 
and a-enhancement, for galaxies in low- and high-density 
environments. In contrast to the works listed above, they 
concluded that the luminosity-weighted age of the bulk of 
ETGs is independent of the environment. Low- and high- 
density environments actually differ because of the fraction 
of galaxies showing signs of ongoing star formation. The 
percentage of "rejuvenated" ETGs is found to increase sig- 
nificantly in the lower density environments. It is not clear 
if this scenario is also able to explain the environmental 
dependence of the FP. As noticed in Sec. 17.11 the varia- 
tion of FP intercept in r band remains essentially the same 
when restricting the analysis to the subsample of ETGs with 
lower contamination from objects with non-genuine ETG 
morphology. Moreover, the FP offset is computed from the 
median values of the distributions of effective parameters 
and velocity dispersions, hence reflecting more closely the 
bulk properties of the ETG's population. 



Another important finding of the present work is the 
lack of any evidence for ETGs in the field to be less metal- 
rich than t hose in the clu s ter en vironment, as found, for 
instance, bv lGallazzi et afl f2006f). The optical+NIR varia- 
tion of the FP offset implies a difference in metallicity of 
—0.043 ± 0.023 per decade in local galaxy density, corre- 
sponding to a difference of ~ 16 ± 8% between ETGs in 
low- and high-density environments. Hence, our FP anal- 
ysis is consistent with field galaxies being more metal-rich 
than ETGs in groups. Rose et al. (1994) found a similar re- 
sult and interpreted this difference as reflecting the fact that 
star formation and chemical enrichment in ETGs in clusters 
were truncated at an early epoch, so ETGs were influenced 
by the environm ent at an early phase of galaxy evolution. 
For comparison, iBernardi et al] (|2006h found no detectable 
difference in the m etallicity of fiel d and c l uster ETGs, while 
other s tudies, e.g. IThomas et all l|2005l ); Ide La Rosa et aD 
l|2007h :[ci emens et al] (|2009l ) found evidence for galaxies in 
dense environments to be less metal rich than those in the 
field. 

Different from most of the works listed above, we do not 
only analyse the average difference between the properties of 
group and field galaxies, but also follow the variation of FP 
coefficients as a function of local galaxy density. We find that 
the offset of the FP smoothly changes from the highest den- 
sity regions, in the cores of parent galaxy groups, through 
the outskirt group regions and the field. Such variation is 
caused by a variation in age, i.e. ETGs have progressively 
younger stellar populations from the cluster core s through 
the fie ld, in agreement with the recent results of iBernardil 
(2009) and PGF10, who found that, for a given parent halo 
mass, central galaxies are systematically older than satel- 
lite galaxies. In contrast to PGF10, we find a positive (or 
eventually no) gradient in galaxy metallicity as a function 
of local density, that would imply centrals to be either less 
metal-rich or as metal-rich as satellites (see above). 

Our findings can be compared to the expecta- 
tions of_semi^aiiah£tw_models (SAMs) of galaxy forma- 
tion (|De Lucia et alj|2006h (hereafter deL06). The model of 
deL06 predicts an age variation, from cluster cores to the 
field, of ~ 2 Gyr (see their fig. 8). The trend is mostly due 
to the fact that, in a hierarchical scenario, halos in a re- 
gion of the Universe doomed to become a cluster are those 
collapsing earlier and merging faster. This would also imply 
that the age variation as a function of local density might be 
larger for more massive halos, in agreement with our finding 
that the variation of the FP offset with the local environment 
is larger for galaxies residing in clusters, relative to groups. 
Over the large environmental range spanned by the SPI- 
DER sample (three decades in local galaxy density) , we find 
an age variation of S log* = 0.14 ± 0.02, i.e. St/t ~ 32 ± 5% 
(~ 3.2±0.5 Gyr, for a galaxy formation epoch of 10 Gyr). For 
a pure age model of the variation of FP offset with environ- 
ment (see Sec.JS]), this difference reduces to 8t/t ~ 28 ± 3% 
(~ 2.8 ± 0.3 Gyr). Both values of St/t are larger than that 
predicted from the SAM, but not much different considering 
all assumptions and different definitions of environments in 
both cases. We can also notice that SAMs predict field galax- 
ies to be less-metal rich than their counterparts in groups, 
in contrast to our FP analysis. As discussed by PGF10, cur- 
rent SAMs still do not provide an accurate description of 
galaxy's metallicities: several recipes, such as the treatment 
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of SN feedback (see e.g. iBertone et alj|2007f l. tidal interac- 
tions, and recycling of SN ejecta, might be able to improve 
the comparison with the the environmental trends seen in 
the data. 



11 SUMMARY 

We have analysed the environmental dependence of the 
FP relation of ETGs, combining optical (SDSS) and NIR 
(UKIDSS) data. Environment is defined using the largest 
group catalogue, based on 3D data, generated from SDSS 
in the low-redshift-Universe (0.05 ^ z ^ 0.095). The main 
results can be summarised as follows. 

• The intercept, "c", of the FP significantly changes with 
environment. Local galaxy density is the main driver of the 
environmental dependence, with "c" decreasing from low to 
high density regions. Weaker correlation is observed between 
"c" and the mass and mass density of parent galaxy groups. 
These latter trends are entirely explained by the change in 
the average local density of galaxies residing in different bins 
of (halo) mass density. The behaviour of "c" with local den- 
sity becomes steeper as the parent group mass increases. 

• The variation of "c" with local environment implies that 
galaxies in higher density regions have on average a larger 
mass-to-light ratio. Since the typical dynamical mass does 
not change significantly with environment and the trend of 
"c" is similar in all wavebands, from g through K, we con- 
clude that the difference in mass-to-light ratio is mainly due 
to a difference of galaxy stellar population properties. 

• The variation of "c" implies that the luminosity- 
weighted age smoothly increases with local density by about 
0.048 dex (~ 11%) per decade in local density, while metal- 
licity tends to decrease. Field galaxies are, on average, 18% 
younger than their cluster counterparts. As galaxies in our 
sample span a range of three decades in local density, the 
age difference between galaxies in the highest density clus- 
ter regions and the field amounts to ~ 32%. We have also 
found (at 2.5 a) that the variation in age per decade of local 
density might be larger, up to a factor of two, for galaxies 
residing in more massive groups. 

• The slope "a" is smaller (at ~ 3 a) for galaxies in groups 
compared to the field ones. Consistent with that, we also find 
a weak (2.5 a) correlation in the sense that a decreases as the 
local density increases. "6" tends to become higher when Ejv 
increases, which is mainly due to galaxies inhabiting more 
massive clusters, as for poor groups no trend is detected. 
The difference in "a" between group and field galaxies seems 
to be similar in all wavebands, from g through K, maybe 
implying that the different tilt of the FP between group and 
field galaxies is related to a different variation of dark-matter 
fraction and/or non-homology along the mass sequence of 
ETGs. 

• The waveband variation of "a" is very similar in all the 
environments: it increases, by ~ 16%, from g through K. "b" 
varies by less than ~ 3% with waveband. For field galaxies, 
"6" increases by ~ 2.5% from g through K (at 3 a), with the 
difference in "6" between group and field galaxies becoming 
smaller in the NIR. For galaxies in more massive parent ha- 
los, no variation of "6" with passband is detected at all. 
These findings may imply that the correlation between stel- 



lar population properties and mass is somewhat shallower 
for galaxies residing in higher density environments. 

As a final remark, we can look at the FP of ETGs as a tool to 
probe physical mechanisms shaping galaxies as we observe 
today. The environmental variation of its coefficients may be 
interpreted in the following way: "a" - reflects some effect 
like homology or dark matter content; "b" is connected to 
how the stellar populations vary with mass; and "c" shows 
how the bulk of system's luminosity varies. These rules can 
be applied in the study of ETGs as a function of environ- 
ment, and at higher redshift, helping us to comprehend, even 
if just in a phenomenological fashion, how galaxies form and 
evolve. 
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